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ABSTRACT 

Context. The rapid neutron-capture process, which created about half of the heaviest elements in the solar system, is believed to have 
been unique. Many recent studies have shown that this uniqueness is not true for the formation of lighter elements, in particular those 
in the atomic number range 38 < Z < 48. Among these, palladium (Pd) and especially silver (Ag) are expected to be key indicators 
of a possible second r-process, but until recently they have been studied only in a few stars. We therefore target Pd and Ag in a large 
sample of stars and compare these abundances to those of Sr, Y, Zr, Ba, and Eu produced by the slow (s-) and rapid (r-) neutron-capture 
processes. Hereby we investigate the nature of the formation process of Ag and Pd. 

Aims. We study the abundances of seven elements (Sr, Y, Zr, Pd, Ag, Ba, and Eu) to gain insight into the formation process of the 
elements and explore in depth the nature of the second r-process. 

Methods. By adopting a homogeneous one-dimensional local thermodynamic equilibrium (ID LTE) analysis of 71 stars, we derive 
stellar abundances using the spectral synthesis code MOOG, and the MARCS model atmospheres. We calculate abundance ratio 
trends and compare the derived abundances to site-dependent yield predictions (low-mass O-Ne-Mg core-collapse supernovae and 
parametrised high-entropy winds), to extract characteristics of the second r-process. 

Results. The seven elements are tracers of different (neutron-capture) processes, which in turn allows us to constrain the formation 
process(es) of Pd and Ag. The abundance ratios of the heavy elements are found to be correlated and anti-correlated. These trends 
lead to clear indications that a second/weak r-process, is responsible for the formation of Pd and Ag. On the basis of the comparison 
to the model predictions, we find that the conditions under which this process takes place differ from those for the main r-process in 
needing lower neutron number densities, lower neutron-to-seed ratios, and lower entropies, and/or higher electron abundances. 
Conclusions. Our analysis confirms that Pd and Ag form via a rapid neutron-capture process that differs from the main r-process, 
the main and weak s-processes, and charged particle freeze-outs. We find that this process is efficiently working down to the lowest 
metallicity sampled by our analysis ([Fe/H] = -3.3). Our results may indicate that a combination of these explosive sites is needed to 
explain the variety in the observationally derived abundance patterns. 

Key words. Stars: Population II - Stars: abundances - Supernovae: general - Galaxy: Halo, chemical evolution 



1. Introduction 

The heavy elements beyond the iron-peak are not created in the 
same way as the lighter elements, many of which form via hy- 
drostatic core or shell burning in the star. These elements are 
generally created by various neutron-capture processes taking 
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place as either the result of mixing in very evolved stars or ex- 
plosions^. 

Previous studies have shown that the slow neutron-capture 
(s-) process can be classified into two sub-processes, namely 
a weak s-process creating the lighter of the s-process iso topes 
dPrantzos et al.lll990tlHeil et al]l2009r.lPignatari et al.l2010i) . and 
a mai n s-process creating heavy isotopes, s uch as those of bar- 
ium dKappeler et alJl989tlBusso et alJl^lGalhno et alj2006t 
ISneden et al.l l2008). The sites of the rapid neutron-capture (r-) 



1 We disregard proton processes here. 
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processes remain unclear, and the exact conditions under which 
they operate continu e to be investigated. Since the time of 
Burbidge et al.l (1 19571) . it has been evident that an explosive en- 
vironment is needed to provide the proper conditions for an r- 
process to happen. After several attempts to make site-dependent 
predictions of the neutron-capture processes, iKratz et all d 1993b 
provided a site-independent approach using the so-called wait- 
ing point approximation, which is based on the best available 
nuclear physics to shed light on the r-process. Nevertheless, 
the conditions are still poorly constrained. A number of sites 
have b e en suggested: neutron st ar mergers (tFre iburghaus et alj 
Il999bt iGorielv et all 1201 lallbk IWanaio & Jankal 120121) mas^ 
sive core-collapse supernovae (SNe) dWasserburg & OianCOOOt 
Arga st et all 120 04). neutr i no-driven winds from core-collapse 



19941 



SNe dDuncan et al.l I l986t iMeveri 119931: [Takahashi et al.l . 
Wooslev et al. 1994; i Freiburghaus et al.l Il999at IWanaio et al.1 
200UlFarouqi et al.ll2009ll2010HArcones & Montesll201 lb. low- 
mass SNe fr om collapsing O-Ne-M g cores (Wanaio et al. 2003; 
or iron cores ISumivoshi et al.ll200lT) . However, no consensus on 
the formation site has been reached. 

Observationally, the discovery of r-process-rich stars which 
contain a factor of 20-100 more heavy elements than nor- 



mal Population II halo stars; see Hill et al. 


2002:ISneden etal. 


2003 




Christlieb et al. 2004; Barklem et al. 


120051 iFrebel et al. 


2007 


Havek et al. 2009; Aoki etaD 2010t |Cowan et alj 201 1 



has offered an important opportunity to study in greater de- 
tail the r-process and its characteristics. By comparing light to 
heavy neutron-capture e lements (i.e. 38 < Z < 50 vs Z > 56) , 
some of these studies (Snedenetal. 2000; Westin et al. 2000; 



Johnson & Bolte 2002; Christlieb et al. 2004; Honda et al. 2004 
B arkl em et al l 120051: Honda et alj 120061 120071: IFrancois et al 

" l2008bl: IRoederer et al 



2007[ ISneden et al.l 120081: iKratzetal 
2010) have revealed a departure of the "light' 



neutron cap- 
ture elements from the main solar-scaled r-process distribution 
curve, which was interpreted as an indication of an extra pro- 
cess. This suggests that the r-process may als o split into two sub- 
channels, namely a 'weak' and main one dCowan etal] [19911 
IWanaio & Ishimarul2006HOtt & Krat J2 008). which are respon- 
sible for the production of the lighter and heavier r-process iso- 
topes, r espectively. The nom enclature is used to match the s- 
process (Kappeler et1illl989l) . 

The 'weak' r-process has received a lot of recent attention, 
but is still poorly constrained despite the many attempts that 
have been made to understand this process. Some of the pro- 
posed processes are the lighter element primary p rocess (LEPP. 



iTravaglio et al. 20041; lArcones & Montes) 
process (IKratz et~ T 20071: iMontes et al. 



2011) the weak r- 
2007; Farouaietal. 



2009; Wanaio et al.ll201 lb . the vp-process (iFrohlich et al.l l2006). 



and sever al more processes a nd comparisons, which ca n be 



found in ICowan et alJ l200ll lOian & Wasserburgl 1200 ll and 



ISneden et al.l 120031 These processes can be considered when 
attempting to explain the abundances of the lighter heavy el- 
ements, which have been found to deviate from the solar- 
scaled r-process patterrfl Palladium and silver are among these 
light er heavy elements. S ilver was studied for the first time 
bv ICrawford et al l (Il998h more than a decade ago in a small 
sample of metal-poor stars. They applied a different hyperfine 
split oscillator strength from the one we adopt here, which 
together with the higher solar Ag abundance helps us to ex- 
plain the low silver abu ndances they derive. A few years later, 
Uohnson & BoTtel d2002) studied both Pd and Ag in a sample 
that is the only other relatively large sample where both Pd 



and Ag were analysed. H ence, we compare our results to theirs. 
lHansen & Primasl (1201 lb presented the first results of an anal- 
ysis of Ag and Pd in a large sample (55 stars) that demon- 
strated the need for an extra production channel. Here, we ex- 
tend the study to the entire sample (71 stars) and compare our 
derived Ag and Pd abundances to those of five other heavy el- 
ements, namely Sr, Y, Zr, Ba, and Eu. We furthermore wish 
to explore the nature of the second r-process in depth by in- 
vestigating the trends of two particular tracer elements, palla- 
dium and silver. We characterise and constrain the fundamen- 
tal parameters of the formation of these elements by means of 
a detailed comparison to yield predictions from several of the 
above-mentioned astrophysical sites and objects. Silver and pal- 
ladium are important for two reasons. First, silver is predicted 
to be a good tracer of the weak r-proces s since nearly 80% of 
its solar system abundance is predicte d (lArlandini et al.l Tl999: 
ISneden etali20 08: Lod ders et al.l2 009) to have come from the r- 
process, and more than 7 1 % of the r-process is estimated to orig- 
inate from the weak r-proc ess (IKratz et al.ll2008bt iFar ouai et al. 
2009: IRoederer et alJuOlol) . For comparison, only 54% of palla- 
dium is created by the r-process (lArlandini et al.lll999b . Second, 
these two elements h ad only been stud i ed in a small number 
of stars (< 20) until Hanse n & Primasl d201 lh . whereas many 
other neutron-capture e lements such as Ba have been studied in 
hundreds of stars (e . g. iReddv et al.ll20"06t iBarklem et aD 120051: 
IFrancois et al.ll2007t lRoedererll2009b . A study of palladium and 
silver provides astrophysical information on a poorly studied 
part of the periodic table. 

The paper is organised as follows: Sect. 2 describes the ob- 
servations and data, Sect. 3 outlines the stellar parameter deter- 
mination, Sect. 4 presents new atomic data and calibration of 
the line list, Sect. 5 presents the abundance analysis, and Sect. 6 
and 7 provide the results and discussions of our abundance and 
model comparisons, respectively. Finally, our conclusions can be 
found in Sect. 8. 



2. Observations and data reduction 

Our sample consists of a mixture of dwarf and giant stars, which 
were observed at high resolution (R > 40,000). The dwarfs 
were observed in the years 2000 - 2002 with the Ultra Violet 
Echelle Spec t rograp h at the Very Large Telescope, UVES/VLT, 
iDekker et al.l d2000h for a project targeting beryllium, which 
requires high signal-to-noise (S/N) data of the near-ultraviole t 
(near-UV) particularly the Be doublet at 313 nm (Primas 2010). 
Similarly high quality data are also needed to detect silver and 
palladium (328.6, 338.3 nm, and 340.4 nm, respectively). The 
spectra cover the wavelength ranges ~ 305 - 680 nm (in some 
cases up to 1000 nm), including the wavelength gaps between 
the CCD detectors. All of our UVES spectra have a S/N > 
100 per pixel at 320 nm. The dwarf spectra were reduced with 
the UVES pipeline (v. 4.3.0). The pipeline performs a standard 
echelle spectrum data reduction. It starts with bias subtraction, 
removes bad pixels due to e.g. cosmic ray hits, and locates the or- 
ders. Then a background subtraction is followed by flat field di- 
vision, order extraction, and wavelength calibration, and finally 
the orders are merged. We tested the quality of the dataproducts 
against a manual data reduction carried out in IRAF|J because 
previous versions of the pipelines had problems with the order 



Solar-scaled r-process abundance: N r = N - N s 



3 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
of Research in Astronomy, Inc., under contract with the National 
Science Foundation. 
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merging. However, this pipeline performs very well and the re- 
duced data were compatible with manually reduced data. Finally, 
the reduced spectra were radial velocity corrected/shifted via 
cross correlation, coadded, and had their continua normalised 
(in IRAF). 

The spectra of the giants were instead extracted from pub- 
lic data archives of the Very Large Telescope (VLT) and the 
Keck telescopes. In both cases, the spectra were observed with 
the hig h-resolution spectro graphs available on both sites, i.e. 
UVES dDekker et all l2000h on th e VLT and High REsolution 
Spectrometer HIRES dVogt et al.l 11994 on Keck. The wave- 
length coverage of HIRES spans 300 - lOOOnm, which is very 
similar to the wavelength range of UVES but might have gaps 
above 620nm. Only spectra of high and comparable (to the 
dwarfs') quality were added to the sample. The giant spectra ex- 
tracted from the respective archives had already been reduced, 
and were carefully inspected, radial velocity shifted, coadded, 
and continua normalised in IRAF. 



Sample 

The final stellar sample consists of 42 dwarf and 29 giant field 
stars, belonging to the Galactic halo, the thick, and the thin disks. 
The sample spans a broad parameter range exceeding 2000 K in 
temperature, 4 dex in gravity, and 2.5 dex in metallicity. Such a 
sample composition allows us to explore the chemical evolution 
of the Galaxy, as well as test the different chemical signatures of 
different stellar evolutionary stages. This in turn can shed light 
on the importance of mixing and non-local thermodynamic equi- 
librium (NLTE) effects. 

Our sample includes some of the most well -known r-proces s 
enhanced giant stars including CS 310 82-001 dHill et al.ll2002l) . 
which we comp are to CS 22892-0 52 dSneden et al.ll2003l) . and 
BD +17 3248 dCowan et all 1201 lb . We note that only one r- 
process enhance d metal-poor dw arf star has been found and 
observed so far dAoki et al.l 1201 Oh . which is not included here. 
Furthermore, silver lines can be detected in giants of all the 
metallicities studied here, but can only be detected in dwarfs 
with [Fe/H]>~ -2.0. This may introduce a small sample 
bias towards metal-poor r-process enhanced giants. No carbon- 
enhanced stars were included in our sample. 



3. Stellar parameters 

We followed different methods to determine the optimal set of 
stellar parameters. With such a large sample, we faced some 
difficulties in applying the same method to the determination 
of the stellar parameters for the entire data-set. The effective 
temperature of most of our stars was derived from colour- 
T e g calibrations to which we applied the necessary band-filter 
and colour corrections. In thi s respect, we tested s everal differ- 
ent co lo ur calibrations from Alonso et al.l (fl996l). lAlonso et all 
(19991). iRamirez & Melendezl d2005l). iMasana et all d2006), 
Onehagetal. (2009i). and Casagrande et alT i 20101) . who make 
use of both ( V-K ) and (b-y) colour indices. In the end, we chose 
the calibration of lAlonso et al.l d 19961, [1999) because these lead 
to temperature predictions that generally fall in the middle of 
the range shown in Fig.Q] The temperature has a large influence 
on the derived stellar abundances. Hence, we wished to avoid 
systematic effects in the abundances by over-/under-estimating 
the temperature, and therefore selected an intermediate temper- 
ature scale. The photometry was from 2MASS (K) and Johnson 
V (the V-K was taken from lCutri et al.ll2003l) and the parallax 



was taken from the Hipparcos catalogue dPerrvman et al.llT997h . 
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Fig. 1. Effective temperatures derived for eight stars (of differ- 
ent metallicity, from higher to lower as one moves from left to 
right along the x-axis) with seven different colour- T e s calibra- 
tions (see figure legend). 



Our final effective temperatures are based only on the (V-K) 
colour index. Among the indic es, we considered i t to be the most 
metallicity-independent one dAlonso etalJ 119991) . since infra- 
red magnitudes are less affected by reddening (K is the only 
infra-red magnitude that is available for all our sample stars). 
Additionally, the temperatures derived for the dwarfs based on 
this colour are in good agreem ent with those determined via H^ 
line fitting dNissen et al.l 2007). We note, however, that the (b-y) 
colour tends to predict slightly higher temperature values than 
(V-K). 

The reddening corrections, E(B - V), were mos t ly de- 
rived from the Schlegel dust mapfl (Schle gel etalJ |l998) 
and corrected according to iBonifacio et al.l (l2000l) if they 
exceeded 0.1 mag. For a few stars, w e took the corresr. 



E(B - V) values from the literature (Niss en et al.ll2002l [2004 
120071) . We applied the formula of lAlonso et al.l (Il996l) of E(V- 
K ) _ 2.12E(B-V), which corr e sponds to the average of those 
oflRamfrez & M elendez d2005l) . lKinman & Castellil d2002l). and 
Nisse n et al.l (l2002h . A filter conversion of -0.04 from |Bessell 
(2005, 2MASS to Johnson) transformed the K magnitudes 
from the 2MASS to the Johnson system, and brought both 
magnitudes to the Johnson scale leading to: 



V-K 



OJohnson 



Vjohnson - ^2MASS - 0.04 - 2.12E(B - V). 



Having both magnitudes on the Johnson scale, we converted 
V-K from Johnson to TCS (Observatorio del Teid e), which can 
be do ne by applying the following relation from lAlonso et al.l 
d!994h 

(V - K) TCS = 0.05 + 0.994(V - K) ]ohnson . 

This last part of the filter conversion - Johnson to TCS - cor- 
responds on average to +0.04 mag. We keep all transformations 
for the sake of accuracy. 



http://spider.ipac.caltech.edu/starr/jarrett/irsa/dust.html 
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In the case of stars (typically, from the disk) found to have 
unrealistically large E(B - V) values we decided to derive their 
temperatures spectroscopically. The gravity was calculated from 
Hipparcos parallaxes by applying the classical formula 

log A = log £ + 4 log ^ - 4V + 0.4BC + 2 log n + 0. 12, 

where M is the mass, Vo is the dereddened apparent magnitude, 
BC is the bolometric correctiorQ, and n is the parallax. Stellar 
masse s were taken from the literature (iNissen et al.ll2002l [2004, 
120071) . On the basis of lAlonso etalJ (Il995h . we calculated the 
BC for each of our stars. For the few stars for which no parallax 
was available, we constrained their gravities by enforcing ioni- 
sation equilibrium between Fe I and Fe lfj In general, the Fe I 
and II abundances are in good agreement, although we note that 
for stars labelled with an 'a' or 'c' in online Table IB1 the par- 
allax was neglected (either owing to their large uncertainties or 
wide-ranging Fe I and II abundances when the gravity was de- 
rived from the parallax). The metallicity was derived from Fe I 
equivalent widths (EWs) and is in good agreement with previous 
studies. The microturbulence was determined by requiring that 
all Fe I lines yield the same abundance regardless of their EW. 
The final values and adopted methods are presented in the online 
material [B] 

3.1. Error estimates for the stellar parameters 

The largest source of uncertainty in estimating the temperature 
is the dereddening of the colour indices, e.g. applying overes- 
timated reddening values from the Schlegel dust maps to stars 
close to the Galactic plane. These can easily translate into un- 
certainties of several hundred Kelvin in the derived tempera- 
ture. Disregarding these extreme cases, we found that the gen- 
eral uncertainty in the reddening values is usually 0.05 mag. 
Combining these values of 0.05 mag with the uncertainty due to 
the Johnson-2MASS transformation led to typical uncertainties 
of the order of 100 - 150 K. A slight magnitude-temperature off- 
set was found between giants and dwarfs owing to the stronger 
colour dependence of the dwarfs' temperature compared to that 
of the giants. Similar uncertainties were found for the excitation 
temperatures. 

Since all stellar parameters to some extent are inter- 
dependent, we also tested the influence of gravity and metallicity 
on the temperature. For instance, an uncertainty of +0.15 dex in 
metallicity has a negligible effect on the temperature (the uncer- 
tainty is usually a few Kelvin). An uncertainty of ± 0.2 dex in 
gravity causes an uncertainty in the temperature of < +1 - 10 K. 
Finally, the microturbulent velocity is found to have a negligible 
impact on the temperature. 

The main uncertainty in the gravity comes from the uncer- 
tainty in the parallax, which is on average +1.0" (iPerrvman et al.l 
Il997i) . This translates into < 0.2 dex in logg . A change of +100 
K in temperature only causes a gravity uncertainty of +0.04 dex. 
By altering the gravity by ± 0.25 dex, the Fe II abundance is 
affected by ±0.15 dex, whereas the Fe I abundance remains ba- 
sically the same. 



5 Adopted from lNissen et al.l d!997h 

6 In total, we have 13 stars for which no reliable information on ei- 
ther their (V - K) colour, parallax, or reddening correction, E(B - V), is 
available. Hence, we resort to spectroscopically derived stellar parame- 
ters, i.e. excitation temperatures and gravities constrained via Fe I/Fe II 
ionisation equilibrium (see also letter 'a' and V in online material B). 



The metallicity is based on EW measurements for which Fe 
I and Fe II lines provided consistent results, usually agreeing 
to within 0. 1 dex. Since our derived metallicities closely match 
those found in the literature (most of our stars are well-studied 
Galactic halo and disk stars), our typical adopted uncertainty in 
the metallicity is ±0.1 dex (±0.15 dex in only a few cases). 

For the microturbulence velocity, we estimated uncertainties 
of the order of 0.15 km/s, stemming from the uncertainty in the 
[Fe/H] and the uncertainty in the Fe EW measurements (which 
is of the order of ±2 mA, as tested via repeated independent 
measurements). 

4. Atomic data and line lists 

This section is divided into two. The first part presents the newly 
calculated log gf values of silver, and the second part describes 
the adjustments and calibrations carried out on the line lists. We 
first note that similar calculations are not necessary for palla- 
dium. This element has six naturally occurring stable isotopes 
(102, 104, 105, 106, 108, 110), of which only four are accessi- 
ble to the r-process. 105 Pd is the only odd-mass isotope with nu- 
clear spin 5/2 for which hyperfine splitting exists. The effect on 
the oscillator strength is, however, minor, since this isotope only 
contributes 22.33 < %Qof its solar elemental abundance. Hence, we 
continue focusing only on the hyperfine structure (hfs) of silver. 

4.1. Atomic data 

This section focuses on the derivation of the hfs of the resonance 
lines in Ag I. 

Silver has two stable isotopes with mass numbers 107 and 
109, respectively. The nuclear spin is 7=1/2 for each of the iso- 
topes. As a consequence, each fine structure level is split into 
two hyperfine levels. The resonance lines in Ag I connect the 
lower 5s level to the 5p levels. 

The isotopic and hyperfine structures commonly used in 
abundance studies o f the Ag resonance lines are those given in 
iRoss & Allerl (|1972|) . They derived log gf values for the differ- 
ent hyperfine and isotopic components using t he experimental 
studie s of t he relative hfs pattern c onducted by Jackson & Kuhnl 
dl937l) and lCrawford et all (fl949). These are intensity measure- 
ments of different componen ts studied by interferometric exper- 
iments. |RossjfcAllei] JT973) label four components, i.e. two hy- 
perfine components for each isotope. The expected number of 
components are three for each of the isotopes 107 and 109 (see 
TablefTland lA.U . The uncertainty in the old intensity measure- 
ments resulted in a misinterpretation and misidentification of the 
components. 

We derive new hyperfine transition components based on 
several experimental measurements of the hfs from more re- 
cent studies, using the theory of the addition of angular mo- 
menta to derive the hyperfine components. We also derive ex- 
perimental oscillator strengths, log gf values, for the different 
components. The transition energies are derived from unresolved 
Fourier transform spectroscopy (FTS) measurements. 

Hyperfine structure components 

The splitting due to the hfs of a level is given by 

A£m s = ^A his [F(F + 1) - J(J + 1) - /(/ + 1)], 
7 http://www.tracesciences.com/pd.htm 
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where Am s is the hyperfine magnetic dipole constant. For nuclei 
with larger spin, the electric quadrupole moment can be signif- 
icant, but for nuclei of spin 1=1/2 , as for Ag, only the magnetic 
dipole is non-zero dCowan|[l9 8li[). The quantum numbers I, J, 
and F are those related to the nuclear spin, total angular mo- 
mentum of the electrons, and the total angular momentum with 
the nuclear spin taken into account, respectively. This expression 
assumes that the coupling among the electrons, resulting in a to- 
tal angular momentum J, is much stronger than the coupling to 
the nuclear angular momentum /. The interaction between / and 
J are coupled to a moment F. 

The energy splitting for a given level can thus be derived 
from the hyperfine constant. The hyperfine c onsta nts Aj,f s for the 
5p levels were measured bv lCarlsson et al.l (1 19901) by observing 
quantum beats. The splitting of the 5s level is an order of mag - 
nitude larger and was measured by Dahmen & Penselin (1967). 
From the energy splittings, the relative wavelengths for the tran- 
sitions can be derived. 

The intensity ratios for the transitions between the different 
hyperfine components can be derived using the exp ressions for 
the addition of angular momenta (e.g. lCowanlll981l) . where the 
decay in each channel is proportional to 



A oc (2F+ l)(2F' + 1) 



./ 

F' 



and the prime is for the lower level. From the hyperfine constants 
of the 5s and 5p levels, the hyperfine pattern with relative inten- 
sities and splitting can be derived. This gives the relative inten- 
sities and positions of the hyperfine components for one isotope, 
but not the relative shift between the isotopes. 

We used the interferometric observations of [Jackson & Kuhnl 
dl937l) to derive the shift between the two isotopes. The resolved 
components in their measurements were, with the aid of the pre- 
dicted hfs for each isotope, used to derive the isotopic shift. We 
used the resolved components (F u - F/: 1-0) to establish the iso- 
topic shift, which are 0.026 cirT 1 and 0.022 cm -1 for the 5s 2 Si/2 
- 5p 2 Pi/2 and 5s 2 Si/2 - 5p 2 P3/2, respectively. The resulting 
structure for the two resonance lines are shown in Fig. |2] 

The absolute wavelengths of the different components were 
derived from the centre of gravity of the resonance lines mea- 
sured bv lPickering & Ziliol (1200 lb . who used a hollow cathode 
discharge and Fourier Transform Spectrometer. The hyperfine 
and isotopic structure are too small to be resolved in the Doppler 
broadened line profiles. 



Transition strengths 

The derivation of the line structure due to isotopic and hyperfine 
structure above give the relative intensities. To use the transi- 
tions for quantitative studies, we need the absolute values, i.e. 
the oscillator strengths (log gf), which can be derived from the 
radiative lifetime of the upper levels. 

The lifetimes for the upper levels of the resonance transi- 
tions, 5p 2 Pi/2,3/2 w ere measured u s ing a laser induced fluores- 
cence technique by ICarlsson et al.l d!990t) . Since there is only 
one decay channel per level, the transition rates (A) are simply 
given by the inverse of the lifetime as A = 1 /r. 

The absolute transition rates can, combined with the relative 
intensities of the hyperfine components for a given fine struc- 
ture transition as discussed above, give the log gf value for the 
individual hyperfine components according to 
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Fig. 2. The hyperfine and isotopic structure of the resonance 
lines of Ag, calculated using a natural isotopic abundance. 



where A is given in nm and g is the statistical weight. These are 
reported in Table IA.1I 

The hyperfine and isotopic structures of Ag is rather small 
and cannot be resolved in the stellar spectrum. The contribu- 
tion from the different isotopes can thus rarely be measured. To 
handle the different isotopes in the stellar spectrum, it is usually 
assumed that the isotope ratio is the same as the natural abun- 
dance: 51.84% for isotope 107 and 48.16% for isotope 109. It is 
convenient to derive the contribution to the Ag absorption lines 
from the different isotopes, normalising to the isotopic ratio. The 
line parameters for a natural abundance mix of isotopes is given 
in TableQ] It should be noted, however, that the true log gf is an 
atomic parameter for each isotope, which is independent of the 
isotopic ratio, and the values in Table Q] are to be used only with 
a fixed isotopic ratio and for a total Ag abundance. For a strict 
treatment of the individual abundances for the two isotopes, the 
values in Table [AT] should be used. 

Figure [3] shows the effect of including hyperfine splitting 
with zero, two or three hfs levels. If we had adopted the log gf 
value available from VALD (the Vienna Atomic Line Databasqfl 
iKupkaFl koOO) without hfs, all the Ag abundances would have 
been overestimated. This effect is even more pronounced in the 
cool metal-rich stars, where the silver lines are stronger. In dwarf 
stars such as the Sun, the new hfs predicted log gf values can 
lead to a difference of <+0.2 dex in the derived silver abun- 
dances, compared to the results based on iRoss & Allerl (1 1972b 
values (see Fig. [3}. Hence, neglecting hfs would lead to overes- 
timated silver abundances. 



Silver isotopes 

Based o n measurements of the visual and near-infrared Ag I and 
II lines (lElbel& Fischer! 19621) . silver is predicted to show a rela- 
tively small isotopic shift, which would barely affect the spectral 



gf = 1.499-10 



U A^gA, 



http://vald.astro.univie.ac.at/~vald/php/vald.php 
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Table 1. Model parameters for the silver resonance lines, assuming an isotopic ratio of 51.84% for isotope 107 and 48.16% for 
isotope 109. 
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2o 
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r 1/2 
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total log gf 


-0.334 


107 


2 Sl/2 


2 P3/2 


0-1 


3280.684 


-0.909 


107 


S 1/2 


P3/2 


1-1 


3280.678 


-1.210 


107 


2 Sl/2 


2 P 3 /2 


1-2 


3280.678 


-0.511 


109 


2 Sl/2 


2 P3/2 


0-1 


3280.686 


-0.941 


109 


Sl/2 


P3/2 


1-1 


3280.679 


-1.242 


109 


2 S[/2 


2 P 3 /2 


1-2 


3280.680 


-0.543 










total log gf 


-0.022 



Notes. The transition strengths (reduced log gf) given are not the true log gf, but were adjusted for the natural isotopic ratio. For studies treating 
the isotopes using individual abundances, the data in Table \AA\ should be used. 




3280.4 3280.6 3280.8 3281.0 
Wavelength [A] 

Fig. 3. The Kitt Peak solar spectrum with spectrum synthesis 
computed with different line lists overplotted: VALD's log gf 
without hfs (dotted blue line); our most recent log gf values 
(dash -dotted red line); a nd the old hfs (only two levels) values 
from lRoss & Allerll 1 972| (R & A, dashed green line). 



line at our spectral resolution. We carried out a test for the near- 
UV lines with natural isotopic abundance (which is ~ 48/52% 
for 109/107 Ag) and compared this to two other test cases with 
ratios of 25/75 % and 1/99% for the 109/107 Ag isotopes, re- 
spectively. The actual change in the synthetic spectrum was less 
than the width of the plotted line. Hence, the change in isotopic 
fraction could be seen in neither our high quality spectra nor the 
high-resolution Kitt Peak spectrum of the Sun. 



4.2. Line list 

We now focus on the silver and palladium lines and their atomic 
data, since these elements are the ones that have been studied 
the least. The line list for the Sr, Y, Zr, Ba, and Eu lines is not 
reported here. They include the most commonly used transitions 



of these elements, and can be found in the online material (Table 
EOT) . 

In general, all atomic data were taken from VALD dKupka F.I 

2000), and we cross-checked excitation potentials and oscilla- 
tor strengths (log gf) against the NISTQ (National Institute of 
Standards and Technology) compilation and recent literature, in 
order to get the most up-to-date line list and best possible abun- 
dances. 

From VALD, we excluded all weak lineal i.e. lines with 
excitation potential higher than 4 eV and log gf values smaller 
than -4 dex. These weak lines have no significant influence 
on the continuum, thus do not affect the derivation of the Ag 
abundances. We note th at the same approach was followed by 
Uohnson & Boltel ((2002), which we adopted to be able to make a 
direct comparison to their (the only other) large available sam- 
ple. 

The silver lines are situated at 3280.7 A and 3382.9 A and 
the palladium line used in this study falls at 3404.58 A . In this 
near-UV region, the molecular lines (OH and especially NH) 
make a significant contribution to the spectrum, and all molecu- 
lar line information was taken from Kurucz's database^ In ad- 
dition, we note that this wavelength region suffers from uniden- 
tified transitions. Therefore, one predicted line from Kurucz - 
the 3382.96 A, Fe I line - was included in our final list in order 
to produce a satisfactory synthetic spectrum. 

For the 3280.7 A line, the red wing is severely affected by 
blends from the Zr II and Fe I transitions. By synthesising the 
region around the blue silver line using the derived metallici- 
ties of the stars, we found that the blending Fe line (3280.76 A ) 
in most cases is overpredicted (red line in Fig. |4}. Because our 
sample covers a large range of stellar parameters, we ran several 
syntheses, for a large number of stars spanning our entire param- 



9 http://physics.nist.gov/PhysRefData/ASD/lines_form.html 

10 By adjusting the VALD 'extract stellar' search the minimum log gf 
around the silver lines found is -3.4 dex, whereas using the VALD 'ex- 
tract all' yields a factor of five more lines reaching minimum log gf 
values of -9.7. This large number of weak lines evidently affects the 
continuum placement. 

1 1 http://kurucz.harvard.edu/molecules.html 



6 



C. J. Hansen et al.: Silver and palladium help unveil the nature of a second r-process 




3280.0 3280.2 3280.4 3280.6 3280.8 3281.0 3281.2 
Wavelength [A] 

Fig. 4. The effect of a wrong log gf of the blending Fe line 
(marked by an arrow) shown for HD 121004. The log gf of Fe I 
transition (red line) over-predicts the Fe line strength, resulting 
in an underestimation of the Ag abundance. The synthesis using 
our adjusted Fe I log gf value is shown in blue. 
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Fig. 5. A spectrum of HD 12 1004 with the results of two spec- 
tral syntheses with different log gf values for Zr; in blue we plot 
the results for -1.5 dex and in red for -1.1 dex. This demon- 
strates that a reduction in this zirconium line's log gf value was 
necessary to obtain a better synthesis and the correct silver abun- 
dances. 



eter space with different log gf values for this line. In the end, 
we constrained the value of its transition probability so that it 
gives a reasonable fit to the entire sample. We thus altered the Fe 
I line logg/ value from -2.231 dex to -2.528 dex. An example 
of this procedure is provided in Fig. |4]for the star HD 121004. 
The value listed in the VALD databa se (-2.231) could be found 
in neither NIST nor the Fe line list of lFuhret al.1(ll988l) . 

Furthermore, we note that with this change we were also 
able to derive consistent solar abundances from both silver 
lines. Both solar spectra, the one observed with and 
the Kurucz Solar Flux Atlas[3 yielded silver abundances that 
differed by ~ 0.3 dex, with the bluer of the two lines giv ing 
the lowest silver abundance. The Kitt Peak solar spectrurro 
which has the highest resolution (R ~ 840,000), also yields dif- 
ferent abundances, of the order of 0.19 dex. The alteration of 
the Fe log gf to -2.528 dex led to an agreement between the 
two Ag lines/abundances within 0.04 dex of the two solar sil- 
ver abundances and yielded a value of 0.93 ± 0.02 dex. This is 
in good agreement w ith the previous solar photospheric abun- 
dances summarised in Asnlu ndet alJ (|2009, where log e(Ag) Q = 
0.94 dex). 

The synthesis of this region requires one more change to pro- 
vide an acceptable fit. Based on equivalent width measurements 
of Zr II lines in the optical (see Sect. |6}, we first determined the 
Zr abundance of each sample star, and used those values when 
synthesising the Ag line at 3280 A. We noticed a similar feature 
as for the above-mentioned Fe line: the Zr abundance derived 
from the Zr II line in the red wing of the Ag line was always 
overestimated by ~ 0.4 dex (in all sample stars) when using the 
Zr abundance derived from the Zr optical lines. We then reduced 
the Zr log gf of the 3280.735 A by 0.4 dex and obtained an over- 
all much better fit (see the blue line in Fig . O . 



12 R ~ 85, 000, http://www.eso.org/observing/dfo/quality/UVES 
/pipeline/solar_spec tram.html 

13 R ~ 500,000, http://kurucz.harvard.edu/sun.html 

14 ftp://nsokp.nso.edu/pub/atlas/fluxatl/ 



There are two additional important blends that contribute to 
the region around 3280.7 A , namely that of Mn I and NH; how- 
ever, for neither of these lines are changes needed to their atomic 
data, but they can be properly synthesised once we determined 
their abundances from other spectral lines/regions. 



1 .0 F 




3382.5 3383.0 3383.5 3384.0 

Wavelength [A] 



Fig. 6. A spectrum of HD 12 1004 (dots) to which two syntheses 
are fitted. The red indicates that the log gf value is too low, while 
the blue shows the properly adjusted log gf for the blending Fe 
I line. 



The 3382.9 A silver line has a strong Fe blend in its red wing 
(3382 985 A). This line is taken from the line list of iMoore etall 
d 19661) . b ecause it was not found in either VALD or NIST. 
However, [Moore et al. ( 1966) only provide the excitation poten- 
tial of this line, and we had to adjust the log gf empirically to 
obtain acceptable fits for this wavelength region. We adopted a 
log gf value of -3.28 +0. 1 dex, which provides a good fit to the 
vast majority of our 71 sample stars. 
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The pall adium line list was part ially based on the line list 
published in lJohnson & Boltd (120021) and partly on VALD. The 
list required few (negligible) empirical adjustments and the so- 
lar value obtained from synthesising the line in the Kitt Peak 
solar spe ctrum was log e ( Pd) = 1.52 dex. As previously 
noted in Han sen & Primasl d201 ll) . this value compares very 
well to the solar ph otospheric abundance of Pd summarised in 
lAsplund et ail ( 120091) . log e (Pd) = 1.57 dex. 

For Ba and Eu, we used the hfs calculated relative o scilla- 
tor strengths from iMcWilliaml (1 1 99 8l) and llvans etail (120061) . re- 
spectively. To derive accurate abundances, we applied a weight- 
ing to the lines from which we synthesised the abundances. For 
barium, we assigned the 5853 A line the highest weight (3) 
since this line is clean, and the 4554 A line has an intermediate 
weight (2) owing to the weak blends. Only when neither of the 
two aforementioned lines were detectable was the 4934 A line 
used (with weight 1 — otherwise it was given a weight 0) ow- 
ing to the severe blends, yielding consistently lower abundances. 
Furthermore, we note that the 4554 A line tends to yield higher 
abundances (~ 0.1 - 0.15dex) than the 5853 A line due to the 
presence of blends. Similarly, we assign weights to the Eu lines: 
4129 A was given the highest weight (3) since it is clean, 4205 A 
an intermediate weight (2) owing to the weak blends, and the 
6645 A line (weight 1 or 0) is only used when the two blue 
lines are neither detectable nor observed. The 4205 A Eu line 
yields abundances that on average 0. 1 dex higher than those of 
the 4129 A line, while the abundances of the 6645 A line agree 
with the 4129 A derived ones. However, the 6645 A line is weak 
and generally only provides upper limits for our stars. 



5. Abundance analysis 

The abund ances were calculated b ased on MARCS model atmo- 
sphere^ (Gustafs son et alJl2008l) . which were interpolated to 
match the stellar parameter s derived for our stars using the code 
written by Mas serorJ(l2006l). Additiona lly, the ID LTE synthetic 
spectrum code MOOG ( Snedenl[l973l version 2009 including 
treatment of scattering) was applied to derive the stellar abun- 
dances. To date, neither NLTE corrections nor three-dimensional 
(3D) model effects have been studied for Ag or Pd. However, 
NLTE corrections can be found in the literature for Sr, Zr, Ba, 
and Eu and we briefly comment on these when we discuss our 
results. 

Owing to the severe line blanketing affecting the near- 
UV/blue part of the spectra of all stars, blending plays a ma- 
jor role, thus spectrum synthesis is required to derive accurate 
abundances of Ag and Pd. Since hfs is substantial for the Ba 
and Eu abundances, we also derived their abundances via spec- 
trum synthesis. For the other elements that we studied (Sr, Y, 
Zr, and Fe), we measured equivalent widths mostly in the red- 
der parts of the spectra to avoid line blends. We measured most 
equivalent widths manually, by fitting Gaussian line profiles in 
IRAF (splot task), e xcept for iron for which we used Fitline 
dFrancois et al.ll2003l) . due to the large number of Fe lines avail- 
able in our spectrio 



15 See http://www.marcs.astro.uu.se/ for model atmospheres in radia- 
tive and convective scheme (MARCS models). 

16 The abundances are calculated as: 

[A/B] = log(A/B) - log(A/B) Q , where log 6(A) = log(^) + 12, 
where N A and N H are the number densities of absorbing atoms of ele- 
ment A and hydrogen, respectively. We adopted a scale where the num- 
ber of H atoms is set to 10 12 . 



5.1. Correlation with stellar parameters? 

To ensure that our abundances are pure tracers of formation and 
evolution processes, and unaffected by spurious analytical ef- 
fects and method biases, it is important to carefully investigate 
the trends of the derived abundances with temperature, gravity, 
and microturbulence. 
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Fig. 7. Abundances of Ag (left) and Pd (right) compared to stel- 
lar parameters. They show a clear division between the dwarfs 
and the giants. No trends could be fitted owing to the very large 

r- 



FigureQshows that no trend with any of the three parameters 
is found, but it is evident that there is an abundance difference 
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between dwarfs and giants. Non-local thermodynamic equilib- 
rium effects could be one possible explanatio n of this difference ; 
other 'possibilities could be mixing effects (ISalaris et al. 1 12000: 
lKornll206"8l : iLind et aD2008l) . microturbulent velocity, an incor- 
rect treatment of the T — r relation in the mode l atmospheres o f 
giants, or unknown line blends in the spectra (lLai et al.l 12008). 
This abundance difference cannot be explained by dif ferences in 
the stellar evolutionary stages (cf. lPreston et al.l l2006). 

The compa rison of the Pd and A g abundances to [Fe/H] 
can be found in lHansen & Primasl d201 lb . where flat trends with 
metallicity were found. This means that the abundances are not 
biased by the stellar parameters or the methods applied to deter- 
mine these, and our abundances can be seen as pure tracers of 
the formation processes. This allows us to apply the abundances 
as direct indicators of the chemical evolution of the Galaxy. 



5.2. Error estimation 

The final error in the derived abundances stems from un- 
certainties in the stellar parameters, the synthesis/equivalent 
width measurements, and the continuum placement. The stel- 
lar parameter uncertainties are (r e jf/log g/[Fe/H]/£): +100K/0.2- 
0.25dex/0.15dex/0.15km/s (cf. Sect. 3.1). Their effect on the 
abundances was constrained by running different models in 
which each parameter was varied by its corresponding uncer- 
tainty, one at a time. 

Furthermore, since we synthesised both Pd and Ag transi- 
tions, we needed to include the uncertainty in the continuum 
placement (about +0.05 dex) and the possible incompleteness of 
stellar model atmospheres, the synthetic code, and the line list 
(i.e. missing atomic data), which all together sums up to an un- 
certainty of +0. 1 dex. Adding all three contributions in quadra- 
ture yields uncertainties of the order of +0.2 dex and ±0.25 dex 
in the Pd and Ag abundances, respectively. The average error in 
the equivalent width measurements of Sr and Y is around 2.5 mA 
and slightly larger for Zr, Ba, and Eu (~ 4 mA). These errors 
were incorporated into the total uncertainty in the abundances 
shown in the figures in Sect. [6] 

Propagating the uncertainties in the heavy element abun- 
dances derived from equivalent width measurements and stellar 
parameters resulted in abundance errors of 0.1 - 0.3 dex. Details 
can be found in TablelC. flandlC.2l 



6. Indications of a second r-process 

To characterise the formation process of Pd and Ag, we com- 
pare their abundances to those of various different elements 
that trace the weak/main s-process and the main r-process. This 
comparison allows us to detect either similarities or differences 
between the yet unidentified formation process of Pd and Ag 
and the known formation processes of the elements we com- 
pare to. For this purpose, we selected the following tracer ele- 
ments, which at solar metallicity are created by the process we 
have listed in Table [2] This means that a correlation of Ag with 
Ba around solar metallicity would indicate that Ag had a com- 
mon formation process to Ba, which in this case would be the 
main s-process. However, at low metallicity this picture changes: 
Sr, Y (and Zr) could be created by char ged particle freeze-outs 
dKratz et alJl2008bt iFarouqi et ai1l2009l) . and Ba mainly by the 
main r-process. We find indications that Zr also receives weak r- 
process con tributions at low ([F e/H] < -2.5) metallicities, which 
agrees with lFarouqi et al.ll2009l ("see also Sect. 17). 



6. 1. Chemical evolution trends of Sr - Eu 

We first compare the elemental abundances of Sr - Eu with FeFl 
to follow the chemical evolution of these elements, and detect the 
onset of the various formation processes. We also compare our 
derived abundances to other studies from the literature, which 
include measurements for some or all of the elements stud- 
ied h ere. The following five large abundance stud i es wer e cho- 
sen: Uohnson & Boltd d2002l J02). IBarklem et all d2005l B05), 
[ Francois et al.l d2007l F07), [Bonifacio ~al.l d2009l B09). and 
koedererl d2009l R 09*). The last (R09*) i s a compilation o f 



previous studies by Edvardsson et al. (1993), Fulbright (200C 
iNissen & Schuster! (1 19971) . and IStephens & Boesgaard (2002 



As mentioned in Sect. [2] we include and compa re with some 
r-proc ess enhanced stars . These are: BD+1 7°3248 dCowan et al.l 
2002) , CS 22892 -052 dSneden et al.ll2003l) . and CS 31082-001 
(Hill et a n i2002l included in our sample). These are clearly la- 
belled in the figures. 



» Giants 
• Dwarfs 

CS 22892-Cfe2 • 




J02 • 
B05 • 

F07 < 
B09 + 



[Fe/H] 

Fig. 8. [Sr/Fe] as a function of [ F e/H| for the entire sam- 
ple, co mpared to [J ohnson & Bolte (2002, J02 — orange as- 
terisk), IBarklem et all d2005l B0 5 — black dots), the 'First 
Stars' gian ts Francois et al. (2007, F07 — green x), and dwarfs 
iBonifacio et alJ d2009l B09 — purple +), respectively. The 
dwarfs from our sample are shown as filled blue circles, while 
filled red triangles represent our giants. Three very enhanced 
stars are shown an d labelled in this and the following fig- 
ures: BD+1 7°3248 dCowan et alj|2002t open black square), CS 
22892-052 dSneden et al.l 120031. filled black triangle), and CS 
31082-001 (IHill et alJl2002i also analysed in this study, hence 
the red triangle). Arrows indicate upper limits to the abundances. 
A flat trend of [Sr/Fe] is seen down to [Fe/H] ~ -2.5, below 
which the scatter becomes dominant. 



Starting with the lightest element, Sr, we see that down to 
[Fe/H] = -2.5, [Sr/Fe] presents a relatively clean and flat trend 
with a mean value around 0.14 dex (see Fig. [8). Below this 
metallicity, the scatter becomes dominant. Only three stars devi- 
ate from this picture (HD175179, HD195633, and G005-040), 
for which only upper limits were attainable from near-UV lines 
(no spectra covering the wavelength range 3800 - 4800 A were 
available in the ESO archive). 

The trend for yttrium is also seen to be flat down to [Fe/H] = 
-2.5 dex (Fig. [9]). We find the same increase in star-to-star scatter 



All abundances are available online — see TablelC. UlC. 2 
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Table 2. We list the elements and the process that they trace at solar metallicity. 



Elements 


Formation process [1] 


Reference 


r-process fraction [8] 


Sr 


85% s-process (weak s-process) 


[1,3,6] 


2.7% 


Y 


92% s-process (in part weak s-process) 


[1,3,6] 


5% 


Zr 


83% s-process (low weak s-process) 


[1,6] 


9.2% 


Pd 


54% r-process (some (<~ 67%) weak r-process) 


[1,2] 


46.9% 


Ag 


80% r-process (mainly (>~ 71%) weak r-process) 


[1,2,5] 


77.9% 


Ba 


81% s-process, (main s-process) 


[1,4,7] 


11.3% 


Eu 


94.2% r-process (main r-process) 


[1] 


94% 



Notes . The references refer to the f ollowing papers ; f 11 [Arlandini etaD d 19991) . [2 1 iFarouqi et al.l J200^). T31 iHeil etal.1 d2009h . r41 lLodders et all 
(2009), r5l lMontes et all d2007l) . r6l lPignatari et al.l (12010). [7] Snede n et al.1 ( I2008I) . and m lBisterzo et alJ ( 1201 ll) for a comparison to more recent 
r-process fractions. 
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Fig. 9. [Y/Fe] vs [Fe/H] similar compa rison samples a s in Fig. [8] 
but also including a fourth sample — Roederer (2009, R09i/o) 
— shown as blue/purple open diamonds indicating stars be- 
longing to the inner/outer halo, respectively. The enhanced stars 
agree with the other comparison samples as well as our sample. 
However, CS 31082-001 is seen to be particularly enhanced in 
Y. [Y/Fe] shows almost no variation with metallicity down to 
[Fe/H] ~ -2.5 dex. 



Zr lines are intrinsically much weaker than, e.g., the Sr and Ba 
resonance lines. 

Figure QT| shows the evolutionary trend of [Ba/Fe] vs [Fe/H] 
which is characterised by a large scatter (> 2 dex) below a metal- 
licity of [Fe/H] = -2.0. The large scatter can be interpreted as 
an indication of different yields from one enrichment event to 
another, creating an inhomogeneous interstellar medium (ISM). 
However, it could also point towards several formation processes 
being at work and releasing very different elemental ratios into 
the ISM. Even whe n correcting the derive d Ba abundances for 
NLTE effects (see lAndrievskv et al.ll2009l) . the scatter is far in 
excess of any possible uncertainty stemming from observations 
and model assumptions. It is therefore a possible indication that 
different formation processes are at play. Figure HTIshows a large 
spread in the europium abundances. 

The evolutionary trends of both [Pd /Fe] and [Ag/Fe] relative 
to [Fe/H] were previously presented in Han sen & Primasl(l201 lb 
and were found to be flat and scattered, similarly to the other five 
elements discussed above. Here, we thus decided to show new 
plots of Pd and Ag abundances, relative to their neighbouring 
elements (see following sub-sections). 

We note that, in general, the r-process enhanced stars fol- 
low the overall trends, but fall on the upper abundance envelope 
as one would expect from their enhancements. For CS 31082- 
001, we re-derived all a bundances and f ound them to agree very 
well with the results of iHill et"a"D (120021) . The only exception is 



of [Y/ Fe] with decreasing [Fe/H] as detected in iRoederer et al.l 
d2010l) . However, the average Y abundance is sub-solar. In gen- 
eral the abundance predictions of the Sr/Y-ratio from SN mod- 
els are found to be very high, most likely due to incorrect solar 
scaled residual^. A too-low solar abundance of Y could have 
explained this, but this does not seem to be the case, since the 
solar phot ospheric and meteor itic Y abundance agree to within 
0.04 dex (lAsplund et alf 2009). making this a trustworthy value. 

The zirconium abundance distribution is also flat and found 
to have a mean value of 0.2 dex down to a metallicity of at 
least -2.5 dex (see Fig.fTOb. The scatter in [Zr/Fe] below [Fe/H] 
= -2.5 is less pronounced than for [Sr/Fe], which may be due to 
there being fewer Zr abundance determinations at low metallic- 
ities compared to, e.g., Sr. One can see from Table IA31 that the 



18 The Sr/Y-ratio can be co rrectly predicted by the high-entropy wind 
models (Farouqi et al. 2009), where these residual assumptions are not 
considered. 
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Fig. 10. [Zr/Fe] as a function of [Fe/H]. Zr does not vary much 
with metallicity. Symbols and colour are the same as in Fig. [8] 
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Fig. 11. [Ba/Fe] plotted vs [Fe/H]. Below [Fe/H] ~ -2.0, a very 
large scatter in all samples is seen. The very large scatter is in- 
dicative of a poorly mixed ISM. Symbols and colour coding as 
inFig.E] 
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Fig. 12. [Eu/Fe] as a function of metallicity. A very large scat- 
ter is seen at all metallicities (also within the various samples). 
Symbols and colour coding as in Fig. [8] 



yttrium, which we propose is caused by uncertainties in the con- 
tinuum placement (+0.1dex) and the profile fitted. The Y lines 
to which we fit Gaussian profiles are very sensitive to the exact 
shape and broadening of the profile, and we can only reproduce 
the observed spectral line by fitting much broader line profiles to 
the Y lines than the surrounding spectral lines. The offset in line 
profile between the Y lines and the nearby other spectral lines 
introduces an 0.3dex abundance offset in our Y abundance. We 
c an attribute our h igher Y abundance compared to that derived 
m lHilletal](l2002h to a combination of uncertainties and offsets. 

The star-to-star abundance scatter revealed by all the elemen- 
tal trends discussed here points to a rather inhomogeneous ISM 
below a metallicity of -2.5 (see Sect. I6.4l for further discussion). 
Below this metallicity, the varying abundances indicate that the 
stars have been affected by different productions (or processes) 
from various nucleosynthetic events. The main contribution at 
these low metallicities must come from primary processes, since 
the sites of the secondary processes (the s-processes) have not 



yet had enough time to both reach the evolutionary stages where 
they yield s-process contributions and have their yields incor- 
porated into later generations of stars. This is why any mon- 
itoring of the r-process is carried out most efficiently below 
[Fe/H] = -2.5. From Fig. [8] — |TT] the s-process might start 
around [Fe/H] = -2.5 dex, since we see a change in the abun- 
dance behaviour (trend flattening/lower scatter) at this metallic- 
ity. Unfortunately, our data do not allow us to identify the metal- 
licity for the onset of the weak s-process, a problem that we dis- 
cuss further in Sect. 16.31 



6.2. Correlations and anti-correlations 

We now turn to a different set of abundance plots, of the type 
[A/B] vs [B/H] (where A and B are two of the neutron-capture 
elements under investigation), to see whether and how they (anti- 
) correlate with each other. This is determined by the abun- 
dance trends to which we fit lines. The slopes determine the 
anti-/correlation. The fitting of linear trends has been made to 
all points (stars) taking their uncertainties into consideration, 
and the uncertainties in the fits are expressed in the figures in 
parentheses. These plots are powerful diagnostics for constrain- 
ing formation processes and can help us to identify similarities 
and differences among the neutron-capture elements. If A and B 
correlate (i.e. the [A/B] ratio is flat across the spanned values of 
[B/H]), it means that they grow in the same way (constant ratio) 
and that they are most likely created by the same process. If they 
anti-correlate (e.g. [A/B] decreases with increasing [B/H]), this 
is usually interpreted in terms of their having different amounts 
of A and B, hence different processes being responsible for their 
formation. To define our terminology, the strengths of the corre- 
lations can be described as follows; a weak/mild anti-correlation 
is stated for slopes between -0.25 and -0.5 and a strong anti- 
correlation is assigned to negative slopes around or steeper than 
-0.5. We choose hydrogen (H) as our reference element because 
we wish to focus only on the formation processes of elements 
A and B. Had we selected iron instead, the interpretation of the 
plots would have become more complex because of the different 
sites contributing to the formation of iron. 

In the following, there are two important factors to bear in 
mind, namely the difference between dwarfs and giants and that 
below [Fe/H] < -2 dex the silver lines could only be detected in 
giant stars. The giants might have been affected by NLTE or mix- 
ing effects, whereas the inclusion of the dwarfs may affect our 
constraints on the formation processes. The giants could be af- 
fected by almost pure r-process yields, whereas the dwarfs might 
carry a mixture of r- and s-process yields. Therefore, we need to 
test the purity of the r-process as we do in Sect. |7] Furthermore, 
it is very important to look for differences in the behaviour of the 
Ag and Pd abundance ratios in dwarf and giant stars (see Sect. 
1641) . 

Now focusing on the formation process of Pd and Ag, we 
start by comparing these two elements to Sr, Y, and Zr, which 
may be formed by the weak s-process elements or charged par- 
ticle freeze-out (depending on metallicity). 

In general, Fig. [13] [T4] and [15] have one common feature, 
i.e. they all clearly show that the elements plotted in each graph 
anti-correlate. Although these anti-correlations are characterised 
by slightly different (negative) slopes, all of these plots agree 
that neither Pd nor Ag are formed by the same mechanism that 
produced Sr, Y, or Zr (i.e. weak s-process or charged particle 
freeze-outs). However, these negative slopes do not merely differ 
randomly between the elements, but there seems to be a clear de- 
creasing trend (i.e. the slopes become shallower) going from Sr 
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Fig. 13. [Ag/Sr](left) and [Pd/Sr] (right) as a function of [Sr/H] is shown here for both dwarfs (filled blue circles) and giants (filled 
red triangles). An anti-correlation is seen in this figure, which is strongest for the dwarfs (see the slopes in the figure). The values 
given in parenthesis are the uncertainties in the linear fits: the first number is the error in the slope, the second number is the 
uncertainty in the intersection with the y-axis. 
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Fig. 14. Left: [Ag/Y] as a function of [Y/H]. Right: [Pd/Y] vs [Y/H]. Legend is described in Fig. [T3l and shown in the figure. 
Anti-correlations between the weak s-process element Y and Ag and Pd are seen in this figure. 




Fig. 15. [Ag/Zr] and [Pd/Zr] vs [Zr/H] to the left and right, respectively. The trend of the fitted line is only slightly negative, which 
could be interpreted as a slight correlation, but the abundances clump. Upper limits to the abundances are indicated by arrows. The 
formulas of the lines fitted are given in the lower left corner for giants and dwarfs, respectively. 
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Fig. 16. An almost fiat trend (correlation) is seen in the figure showing [Ag/Pd] as a function of [Pd/H], which is indicative of a 
similar origin of Ag and Pd. 
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Fig. 17. A strong anti-correlation is seen in this plot of [Ag/Ba] vs [Ba/H] and [Pd/Ba] vs [Ba/H]. Silver and palladium are therefore 
not main s-process elements. 
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Fig. 18. To the left: [Ag/Eu] plotted as a function of [Eu/H], showing a clear and strong anti-correlation. To the right: [Pd/Eu] vs 
[Eu/H]. This means that Ag and Eu are not synthesised by the same process, nor are Pd and Eu. Silver and palladium are therefore 
not produced by the main r-process. 
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to Y and then to Zr. The slopes derived by fitting the data-points 
in [Ag, Pd/Zr] are between -0.37 and -0.18 + 0.07, which thus 
indicate that there is only a mild anti-correlation. We interpret 
this as an indication that Zr may be produced (at least in part) by 
the same formation process producing Pd and Ag. 

When comparing Ag to Pd (see Fig. [T6l. it becomes difficult 
to draw a firm conclusion about the exact trend of their abun- 
dance ratio [ Ag/Pd] as a function of [Pd/H] . Despite the slopes 
overplotted on the graph being indicative of a very mild anti- 
correlation, they may be misleading especially since they take 
into account giants and dwarfs separately. If one were to ignore 
these slopes and consider the entire sample as a whole, we could 
argue that we find a flat [Ag/Pd] trend, especially when con- 
sidering the associated error-bars and excluding upper limits. 
The latter is also supported by our earlier fin ding of an almost 
1:1 li near slope between [Ag/H] vs [Pd/H] ([Hansen & Primasl 
1201 lh . which strongly indicates a common origin for these two 
elements. 

If we now consider how Ag and Pd compare to Ba (Fig. [T7ji, 
which is the most representative tracer of the main s-process, we 
see that both Ag and Pd strongly anti-correlate with Ba, which 
excludes the main s-process as one of the possible production 
channels responsible for the formation of Ag and Pd. At low 
metallicity ([Fe/H] < -2.5 dex), Ba is created by the main r- 
process, which indicates that Pd and Ag are also not created by 
the main r-process, although we compare them to Eu to confirm 
this finding. Finally, Fig. [18] shows that strong anti-correlations 
of Ag and Pd are found with Eu, which means that the process 
forming Pd and Ag cannot be the main r-process. We cannot, 
however, exclude that Ag and Pd are partly produced by the main 
r-process. 

Therefore, the formation process of Pd and Ag is neither 
a charged particle freeze-out, a weak, main s-process, nor a 
main r-process. Both Ag and Pd are seen to form at extremely 
low metallicity ( [Fe/H] < -31 . These results, combined with 
the prediction s of Monte s et alJd2007l) . lKratz et al.1 d2008al) . and 
Farouqi et al. (2009), indicate that their formation process must 
be of primary and likely r-process nature, but we need to re- 
sort to model comparisons in order to characterise this second 
r-process. 

As mentioned at the beginning of this sub-section, one needs 
to keep in mind two caveats when discussing these abundances: 
i) we derived all abundances based on ID LTE model atmo- 
spheres and spectral syntheses; ii) we were able to track the 
evolution of Ag down to the lowest metallicities only with gi- 
ant stars. We adopted the former approach because NLTE cor- 
rections are available for on ly some of the elements inves- 
tigated here, namely S r (e.g. iBelvakova & Mashon kinalll997t 
lAndrievskv et al.1 201 ll and Bergema nn et al, 2012 submitted), 
Zr, Ba (e.g. Andrievsk v et al.l [20091) . and to some extent Eu. 
However, no NLTE corrections have been calculated for our 
two key elements Pd and Ag, and only a few for Y and Zr 
dVelichkoetal]l20Toh . Because we use primarily [A/B] ratios 
(where A can be either Ag or Pd, and B is one of the other 
neutron-capture elements), we decided to keep a ID LTE consis- 
tency across all ratios, instead of correcting only some elements. 
We are, however, fully aware of the importance of NLTE correc - 
tions, and that would ideally be a better way to proceed, were 
NLTE corrections to become available for all elements. As for 
the latter, dwarfs and giants show in general very similar trends 
(see Fig.n"3"1-[T8"I). with the dwarfs having higher abundance val- 
ues than the giants at similar metallicities. However, the over- 
all good agreement between dwarfs and giants suggests that the 



process creating Ag and Pd is likely to be the same at all metal- 
licites. 



6.3. Formation processes and transitions around Zr 

Zirconium and strontium clearly share a common formation pro- 
cess at low metallicities down to and even slightly below [Zr/H] 
= -3 (see the flat correlation for giants in Fig. [T9V A similar 
trend is found when comparing yttrium to zirconium and yttrium 
to strontium. However, at higher [Fe/H] and [Sr/H] abundances 
above - 1 dex, we find an anti-correlation between Sr and Zr for 
the dwarfs. At higher metallicities, this can indicate differences 
in the formation process — or a difference between the process 
primarily responsible for the formation of the two elements. 
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Fig. 19. Zr and Sr correlate in the metal-poor giants which in- 
dicates a similar formation process o f these two elements. This 
is in agreement wi th the findings of iFarouqi et al.l (120091) and 
iKratz et alj (|2008b). At higher metallicities ([Sr/H] > -1) the 
formations of Sr and Zr differ. The upper limits are due to the 
before mentioned lacking visual spectra of the three stars (see 
text). 



Zirconium and barium seem to have different origins, as 
shown in Fig.[20](Zr; e.g. charged particle freeze-out or weak r- 
process vs Ba; main r-proce ss origin at low metalliciti es). These 
finding s confirms those of Faro uqi et al.l 120091 and IKratz et al] 
(2008a, see their Fig. 4), who found a low-entropy charged- 
particle freeze-out process to be the primary formation process 
of Sr, Y, and Zr at low metallicity. Here, we find indications of Zr 
being created in a slightly different way from Sr and Y. Similar 
trends are also seen for [Sr/Ba] and [Y/Ba] ratios, where the gi- 
ants show clear anti-c orrelations. The trend s for giants were al- 
ready reported by e.g. Francois et al. (2007). For the dwarfs, this 
trend is less pronounced and they have a greater scatter in the 
abundances. From the dwarfs' trends, we might conclude that 
around [Ba/H] = -2 the s-process yields from asymptotic giant 
branch stars are no longer negligible formation sites of Ba, and 
that the larger scatter is evidence of multiple formation sources. 
Comparing the giant abundances of Zr to Eu shows that like Pd 
and Ag, Zr is not produced by the main r-process at higher metal- 
licites (see Fig. I2D1 . although we note that Zr and Pd follow a 
weaker anti-correlation with Eu than Ag does. 

In the solar system, Zr appears to have been partly produced 
by the weak and main s-processes (as well as there being a minor 
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Fig. 20. Top: [Zr/Ba] vs [Ba/H] showing anti-correlations. The 
clumping visible at higher [Ba/H] abundances may be indicative 
of some common formation (s-) process for Zr and Ba. Bottom: 
[Zr/Eu] vs [Eu/H], showing a clear, strong anti-correlation over 
the entire range of [Eu/H] values for the dwarfs. This resembles 
the behaviour seen for [ Ag/Eu] vs [Eu/H] . 



contribution from the weak/second r-process), owing to the cor- 
relations (and only mild anti-correlation) of Zr with Sr, Pd, Ag, 
and Ba. At low metallicities, the s-process contribution to Sr, 
Y (and Zr) is substituted with a charged particle freeze-out cre- 
ation. These statements are confirmed in Sect. [7] This means that 
Zr may represent a transition in the periodic table around atomic 
number 40 from the weak s-process/charged particle freeze-out 
process (depending on metallicity) to the weak r-process. This 
second r-process could be responsible for the creation of ele- 
ments in the atomic number range 40 - 50. However, this pro- 
cess would cease to create elements somewhat below barium. 
This upper limit is uncertain owing to the lack of elements inves- 
tigated (observationally in large samples) in the range 48- 55. 
We note that a natural end to the weak r-process from a nuclear 
physics point of view would be around the element tin because 
of the bottle neck occurring at N = 82, beyond which many more 
neutrons are needed to continue the fusion. 



6.4. Discussion 

This section highlights our findings and addresses key points 
mentioned in the previous sections, namely, scatter and inhomo- 
geneities, the presented abundance trends, and differences be- 
tween dwarfs and giants (possibly NLTE effects). 

The consistently large scatter or ISM inhomogeneity seen at 
metallicities below [Fe/H] < -2.5 dex is found in the majority of 
the abundance trends. Many of the large abundance studies have 
found sim ilar large star-to-star scatters at these low metallici- 
ties (e.g.lBarklem et al1l2005t; iPreston et aTll2006t iFrancois etldl 
2007; Bo nifacio etalJl2009h. A NL TE study of the latter car- 
ried out by lAndrievskv et al. (2009), confirmed that the scatter 
in Ba was so large even after applying the NLTE corrections 
to the abundances, that they could not assume that the ISM is 
homogeneous. However, the very low star-to-star scatter of a- 
and iron-peak ele ment abundances provides a coun ter argument 
to this statement dCavrel et al.ll2004t IPreston et al.ll2.006h . since 
these elements would suggest that the ISM is very well mixed. 

Our findings seem to favour an inhomogeneous early ([Fe/H] 
< -2.5) ISM for the reasons that follow. Considering all 
these (alpha, iron-peak, and neutron-capture) abundances above 
[Fe/H] = -2.5, all star-to-star scatters are much smaller and the 
ISM seems to be well-mixed. This implies that single (or a few) 
nucleosynthetic events such as SNe exhibi t smaller effects on 
the st ellar abundances at higher metallicity (Ishimaru & Wanaio 
1999). However, this is not the case below -2.5 dex in metal- 
licity, where we may be witnessing the effects of very dif- 
ferent (single?) explod ing SNe (this was also suggested by 
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pact that various parameters such as peak temperature, mass- 
cut, and entropy have on the SN yields. The a-elements are 
mainly yielded by type II SNe and produced in one process 
only; they do not show this kind of scatter in their abundance 
pattern. The neutron-capture elements, on the contrary, seem to 
have several underlying formation processes, even for the same 
element, which may help explain the variations in the star-to- 
star scatters. The exact site of the neutron-capture elements is 
yet not known, as we have seen in the previous sections, dif- 
ferent ne utron-capture elements might be created via different 
channels (iJohnson & Boltd|2002t IFarouqi et aLll2009l) . Hence, 
the lack of one dominating source could cause a larger scatter 
compared to that of the c-elements. Furthermore, the different 
supernovae that create the neutron-capture elements could, due 
to their differing nature, lead to different neutron-capture pro- 
cesses, i.e. a main and a second r-process, which would help 
us to explain the scatter. Simply put, the inhomogeneity could 
in part be explained by several sources/sites yielding different 
amounts of the neutron-capture elements, whereas the alpha- 
elements are dominated by SNe II which yield relatively similar 
amounts of these ele ments. In contrast to the suggested range o f 
one single r-process dKratz et al. 2008bl iRoederer et all 1201 Ol) . 
which is needed to explain the different abundance patterns of 
HD 122563 and CS22892-052, we confirm that no other group 
of elements be it a, odd-Z, or Fe-peak show this kind of scat- 
ter in abundances when these originate from only one process. 
Furthermore, on the basis of our findings, we see that two r- 
processes (or primary processes) are likely needed to fully ex- 
plain the correlations and the scatter. 

The differences between these processes are clearly evident 
in Fig. [18] where the strong anti-correlation between Ag and 
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Eu, as well as that between Pd and Eu, can be seen. Europium 
is created by the main r-process, a process that requires very 
high neutron number densities to produce Eu (typically around 
1 Q26-28 cm -3 |Kratz et al.ll2007h . whereas the lighter isotopes of 
e.g. Pd can be created in environments with neutron number den- 
sities that are lower by several orders of magnitude. It is impos- 
sible to create Eu in environments with such low neutron densi - 
ties dKratz et alj|2007h iFarouqi et alj|2009l IWanaio et al.ll201 ll) . 
This suggests that the properties of the formation sites for the 
heavy and the light r-process elements differ, or that the pro- 
cesses themselves are different. We note that Fig. [16] indicates 
that the process creates both Ag a nd Pd at almost the same rate 
(see also lHansen & Primasll201 LY The second r-process seems 
to operate effectively at all metallicities down to [Fe/H] = -3.3. 
This process, or its production site, must be less efficient than the 
main r-process. For [Eu/H] > -3, the [Ag/Eu] is below zero and 
rapidly decreases with increasing Eu (see Fig. [TBV However, at 
the lowest metallicities and europium abundances ([Eu/H] < -3) 
the amount of Ag is at the same level or slightly higher than the 
Eu abundance, as can be seen from the [Ag/Eu] abundance being 
larger than zero. This could indicate that the second r-process is 
more efficient at low [Eu/H]. We cannot rule out that Ag and Pd 
both receive small contributions from the main r-process, since 
this process is generally ([Eu/H] > -3) predominant in the ISM 
gas. 

Fi gures [T3l and [T71 show anti-correlations of Ag and Pd com- 
pared to Sr and Ba. At high metallicities, [Fe/H] ~> -1, the s- 
process is far more prevalent in the ISM than the second (weak) 
r-process (e.g. [Ag/Ba] < 0). However, the same figures show 
abundance ratios of around in the range [Fe/H] = -2.5 to —1.0. 
This could indicate that the s-process and the second r-process 
have some features or sites in common (e.g. super AGB stars), 
but this has yet to be confirmed. 

Another important outcome of this study is the discovery 
of Zr as a 'transition' element. Figures [13] to [15] show a grad- 
ual increase in the slopes of Ag and Pd compared to Sr, Y, and 
Zr; i.e. an indication of the growing similarities in their forma- 
tion processes. Within the uncertainties in the slopes, Ag and 
Pd almost correlate/show a very weak anti-correlation with Zr. 
When Ag and Pd are compared to each other (Fig. [TBI , an al- 
most 1:1 correlation is seen. This could be the first observa- 
tional evidence that at higher metallicities ([Fe/H] > -2.5), Sr 
and Y are weak s-process prod ucts, as claimed by iHeil et all 
(120091) and lPignatari et al.l(l2010h (at lower metalli cities Sr and Y 
migh t be created by charged particle freeze-outs; Farouai et al.l 
1201 Ol and IWanaio et alj|201 lb . Zirconium should mainly be an 
s-process element (in the solar inventory), which also receives 
considerable contributions from a type of weak r-process. This 
r-process is responsible for the main production of Pd and Ag. 
The transition from charged-particle freeze-out or weak s- (Sr, 
Y) to 'weak' r-process (Pd, Ag) takes place around Zr (Z = 40), 
hence the name transition element. Moreover, the figures show- 
ing [Ag/Ba] and [Ag/Eu] yield anti-correlations (both strong, see 
Fig.fTTI-fTSTl. suggesting that the formation processes differ. The 
strong anti-correlation with Ba shows that this process is not a 
main s-process and the anti-correlation with Eu demonstrates the 
differences between the main and the second r-process. 

Finally, we turn the discussion to the differences between 
dwarfs and giants. Unfortunately, a full NLTE analysis is not yet 
available, owing to incomplete and complicated model atoms 
of these heavy elements. However, on the basis of previous 
studies of some of the heavy elements such as Sr and B a 
dBelvakova & Mashonkinall997l:lAndrievskv et al.l20TTll2009l) . 
the NLTE corrections can be relatively large for low-gravity 



metal-poor stars. The Sr II abundance may need a correction of 
the order of ~< 0.2 > dex (lAndrievskv et al.ll201 ll) . while the Zr 
II abundance corrections are l ower and gene r ally b etween 0.08 
dex and 0.17 dex according to Velichk o et al.l d2010h . These cor- 
rections are very dependent on the spectral line, the stellar pa- 
rameters, and therefore vary from star to star. Additionally, it is 
insufficient to correct only one of the elemental abundances in 
the abundance ratios we have discussed so far. A detailed NLTE 
study would be needed, but is beyond the scope of this paper. 
Any estimate of the behaviour of the NLTE corrections of e.g. 
silver would be very speculative at this point, although we note 
from Fig. [17] that the [Ag/Ba] ratio of the giants would need an 
NLTE correction of ~ +0.5 dex, as estimated from the offset in 
the figure between the giants and the dwarfs. We note that a frac- 
tion of this estimated value would be due to the NLTE correction 
of Ba. 



7. A comparison to supernova yields 

To gain information on the formation site and process of our 
sample's abundance patterns, we compare these to two differ- 
ent models. The first model (model Ia+b) we focus on is that of 
high-entropy winds (HEW) (iFarouqi et al.l2009[|20Tol) triggered 
by type II SN explosions, whereas the second model (model II) 
is tied to low-mass e lectron-capture SNe (arising from collaps- 
ing O-Ne-Mg cores, IWanaio etalJl201lh . In the last case, the 
neutrino interactions are taken into account. 

To ensure that the abundance-to-model prediction compari- 
son is as informative and complete as possible, we selected eight 
stars distributed in the following way. To probe abundance pat- 
terns that include Ag, two dwarf and six giant stars were singled 
out, where the giant sub-sample includes stars with special pat- 
terns such as r-rich stars. Furthermore, the selection was carried 
out, so that the stars cover a wide range of stellar parameters, 
especially metallicity. By considering stars with a wide range 
of [Ba/Eu] ratios, we include stars with mixed as well as pure 
r-process abundance patterns (see the black diamonds in Fig. 
ED. The stars selected are: CTM153283 and HD106038 (dwarfs 
with mixed r- and s-process patterns), BD +42621 and CS 22890- 
024 (giants; pure r-process tracers), HD 122563 and HD88609 (r- 
poor giants), and CS 31082-001 and HD1 15444 (r-rich giants). 
The dwarf star CD-453283 has a very high europium abundance 
([Eu/Fe] = 0.78), which is the highest Eu abundance measured 
for the dwarf stars in our sample. Over all, this star is overabun- 
dant in elements heavier than Zr. 

IFarouqi et al.l (I20Q9L 120 1 Ot) explored a large parameter space 
especially in entropy, where the values were varied between 
20 and 275 £#/baryon, and the electron fractions, Y e , cover 
the range from 0.4 to 0.49. The wind velocity adopted is 
7500 km/s. The output is neutron-to-seed ratios and correspond- 
ing yi elds/summed abundances. F or further information we re- 
fer to IFarouqi et all (120091 120101) . Owing to the lack of well- 
constrained (3D) supernova explosion parameter output, it re- 
mains unknown whether a high entropy or a low Y e is more likely 
to happen in an actual explosion. Therefore, we carry out two 
different comparisons when contrasting the HEW model predic- 
tions. In model la) Y e is fixed and chosen so that the value repro- 
duces the observationally derived abundances, while the entropy, 
S , is varied. In model lb) the entropy is fixed, while Y e is varied. 
The latter cas e enables a more di rect comparison to the yield 
calculations of Wan aio et al .1 (120111) . 
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Fig. 21. Barium vs Europium. The purity of the process is es- 
timated based on the numbers from Arland ini et al.l (TT999). The 
black diamonds indicate the [Ba/Eu] abundance of the eight stars 
shown in Fig. [24] and [27] 



Model la 



ues of Y e were tested in addition to 0.45. We note that for Y e = 
0.48, the abundances of elements heavier than Zr are underesti- 
mated, whereas Y e = 0.4 predicts too high abundances for these 
elements. The estimates of 0.42 closely reproduce the observed 
abundances, but only the intermediate values of entropies agree 
well with observations — not the high entropies. 

Figure [22] helps us to constrain the entropy value and/or in- 
tervals that provide enough neutron-captures to activate an r- 
process. Our empirical comparison to th e abundances derive d 
for HD 106038 confirms the findings of lFarouqi etall (§009), 
who found two entropy intervals 110<5<150 and 150 < S < 
287 to provide enough neutrons for a weak and a main r-process, 
respectively. Figure[22]shows indeed that the entropies needed to 
create Pd and Ag are likely between 100 and 150 ke/baryon. At 
very high entropies (S ~ 275), no lighter elements (Sr - Ag) 
are efficiently produced, since the fusion continues far past these 
elements owing to the high neutron flux. 

In Fig. [24] we extend the comparison between HEW model 
predictions and derived stellar abundances for eight stars. For 
simplicity, we perform a comparison for only four entropies > 
125^8/baryon. Additionally, we provide in all graphs the [Fe/H] 
and the [Ba/Eu] ratios we derived for each star. The [Ba/Eu] ra- 
tio is especially important, because i t indicates the purity o f the 
r-process (see Fig. [2Tb. According to Arlan dini et al.l (1999). Ba 
is a 81% s-process product, while Eu is a 94% r-process product 
(both in the Sun), hence, the smaller the ratio, the stronger the 
r-process influence is . The r-process i s accep ted as being pure if 



1 111. 1 1 V V IV-l .1 VI LIU. 1 ■ ■ ' IT 111UU\_1.1 Ull. LllUV/j.'J. )J 1 //-■'» J *~~ ' 

Y e , and V exp . All parameters are correlated and define the free the value (- 738) fromlA rndt et al.l201 ll However, stars such as 



The free parameters of the HEW models are entropy, S ~ T 3 /p, [Ba/Eu] < -0_.74 dexlAriandini et al.l (119991). which agrees with 

m A 



neutrons per heavy seed nucleus (F„/F see d)- Neutrino interac- 
tions have only been taken into consideration in terms of their 
estimated impact on the final value of Y e . The model predic- 
tions with a fixed Y e of 0.45 are in good agreement with the 
derived abundances, for both the intermediate (125 £g/baryon) 
and high (200 ke/baryon) entropy (see Fig. [22}. Different val- 



Ye = 0.45 — entropies varied 




Sr 



Zr Pd Ag 
Element 



Ba 



Eu 



HD 122563 dHonda et al]l2006T) provide observationally derived 

upper limits to this range ( 0.2 dex), while a pure s-process 

would be found above 1.14 dex. Furthermore, the metallicity is 
also an important indicator of the predominant formation pro- 
cess, and is therefore included in the figures. Below [Fe/H] = 
-2.5 dex, we generally expect to see r-process yields. In Fig. 
24l - [27] we have scaled all our derived abundances to Sr, since 
we detect this element in most stars and the element is pro- 
duced/predicted at most of the entropies and electron fractions 
considered hereFul We note that plotting [element/Sr] ensures 
that [Sr/Sr] corresponds to zero for all lines. 

Within the error bars, the observationally derived abun- 
dances agree well with the model predictions calculated with 
S = 125 and S = 175fcfi/baryon, although, in a few cases a 
model with S — 150 would have provided closer agreement (see 
Fig. l24l . The neutron-to-seed ratio related to these models are 
in the range Y n /Y see a = 4 - 24^/baryon. From the same fig- 
ure, it is furthermore evident that the heavy elements (Ba, Eu) 
need much higher entropies to be produced. In general, we find 
that the entropy interval facilitating the occurrence of the main 
r-pro cess is 200 - 275 fr g/baryon, which is in good agreement 
wifh[Farouqi et al. (2009). However, we find a slight increase in 
the difference between the weak (125 < S < 175) and the main 
(200 < S < 275) r-process. 

Additionally, these two different processes must be r- 
processes since they are observed in very metal-poor stars and 
show patterns similar to those in the pure r-process stars. 



Fig. 22. A comparison of the abundances of light neutron- 
capture elements as derived in HD 122196 to the yields [X/Sr] 
produced by HEW models. The latter are computed with 
F c =0.45 and varying entropies S (all values listed are given in 
^s/baryon units. The metallicity ([Fe/H]), and [Ba/Eu] ratio for 
the star are shown as well. 



Model lb 

If we now vary the electron fractions, Y e , in the HEW model 
predictions, allowing these to run from S = 2 fc#/baryon to a dif- 
fering final entropy, we see as shown in Fig.[23]that the charged 



Unfortunately, Fe is not predicted in these models. 
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Fig. 23. HEW abundance, Y(Z), predictions as a function of F e for the thr ee different processe s. These predicted abundances are the 
accumulated abundances summed over entropy for a given nucleus (see iFarouqi et alJl2010l for further description). Note that the 
y-axis on the right-side of all three panels indicate the entropy intervals for each of the three processes. Every element is identified 
by a different colour - see legend. 



particle freeze-out, both the weak and main r-processes relate to 
different entropy ranges in the following way: 50 < S < 100 
corresponds to a charged particle freeze-out for a Y e of 0.45, or 
generally speaking, this process takes place when the neutron-to- 
seed ratio is less than one. For the representation of this process, 
we adopted the mid-range value S = 75 fc^/baryon — a value 
that always falls below a neutron-to-seed ratio of one. The weak 
r-process exists at 125 < S < 175 (for Y e = 0.45) or in general 
in the neutron-to-seed ratio from 1 to 15 (where 15 is reached 
at S = 155), the predictions shown in Fig. [23] are made with a 
neutron-to-seed ratio of 5. Finally, the main r-process operates 
when the neutron-to-seed ratio is above 16, and here we have 
shown a ratio of 30 as the representation of the main r-process 
(for a Y e = 0.45, this corresponds to an entropy of ~ 200 — 
similar to what we found in the previous section). The yields in 
percentage for two different electron fractions can be found in 
Table [3] 

In the HEW predictions with Y e = 0.45 and low entropy (S < 
50), mainly iron group elements are produced owing to a very 
low neutron-capture efficiency at these low entropies. Therefore, 
we disregard this part of the entropy range to ensure that we 
produce and consider only heavy elements. Furthermore, not all 
material will necessarily be ejected in the explosion, and some 
fall back is to be expected. 

In the uppermost panel of Fig. |23](the charged particle plot), 
we see that Sr peaks at a Y e of 0.47, i.e. an environment that 
is not very neutron-rich, whereas the Zr yield rapidly increases 
in a more neutron-rich environment and receives contributions 
from both the charged particle process, the weak r-process, and 
a smaller contribution from the main r-process (note the differ- 
ent y-axes in Fig. l23l . This was also seen in Sect. [6] However, 
the contribution from the main r-process was too small to be de- 



Table 3. Percentage of each element created, according to the 
HEW predictions by three different processes. 



Element 


Ch. part. 


weak r 


main 




Y„/Y SCM | < 1 


1 <Y„/Yseed < 15 


Y/i/Y seet | 


Y e = 0.442 








Sr 


97.6 


2.2 


0.1 


Y 


96.7 


3.3 


0.3 


Zr 


81.6 


18.2 


0.2 


Pd 


4.6 


85.3 


10.2 


Ag 


0.7 


82.5 


16.8 


Ba 








100 


Eu 








100 


Y e = 0.493 








Sr 


97.6 


2.3 


0.1 


Y 


94.4 


5.5 


0.1 


Zr 


72.2 


27.6 


0.2 


Pd 


0.4 


72.1 


27.5 


Ag 





63.4 


36.6 


Ba 





1.7 


98.3 


Eu 








100 



Notes. The listed processes are: charged particle freeze-out process (Ch. 
part.), weak r-process (weak r), and main r-process (main r) . These frac- 
tional yields (abundances) are functions of the electron fractions, Y e . 



tected in the abundances. Palladium is according to the HEW 
predictions created by both the weak and the main r-process, but 
as for Zr the contribution from the main r-process is difficult to 
see in the observationally derived abundances (from which we 
found weak anti-correlations between Pd and Eu and Zr and Eu). 
Silver needs considerably more neutron-enriched environments 
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Fig. 24. HEW model yield predictions (computed assuming Y e = 0A5,V expan 
legend) compared to eight different stars (2 dwarfs, in blue; 6 giants, in red). 



7500 km/s and four different entropies, S, see 
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than Zr and Pd, which we see from the slowly increasing 
slope, with decreasing Y e , in the main r-process plot. The heavy 
elements (Ba and Eu) need even more neutrons (lower Y e 's) than 
Pd and Ag to be produced. Comparing Model lb yields to r-rich 
and r-poor stars (CS 31082-001 and HD 122563, respectively), 
we see the increasing need for neutrons with growing atomic 
mass (see Fig. l25l l. In the case of Y e = 0.49, the lighter ele- 
ments can be correctly reproduced by a charged particle freeze- 
out or a weak r-process, although, Ba and Eu require a main 
r-process entropy to be modelled in HD 122563. The environ- 
ment is overall too neutron-poor (or limited to medium entropy) 
to describe the abundances of a r-rich star (see Fig. |25V In the 
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Fig. 25. HEW predictions with a neutron-poor environment 
(Y e - 0.49; upper plot) and a neutron-rich environment (Y e — 
0.442; lower plot) compared to the r-poor star HD122563 and 
r-rich star CS 31082-001. 



Y e = 0.442 case, the lighter to intermediate mass elements are 
within the 0.2-0.25 dex uncertainties correctly reproduced by 
a weak r-process in both stars, whereas Ba and Eu are seen to 
need a main r-process and possibly an even larger neutron-to- 
seed ratio/lower Y e to be correctly reproduced. The need for two 
different processes at work is again expressed by these models 
and r-poor/r-rich stellar abundance patterns. The weak r-process 
cannot create Ba and Eu and the main r-process overproduces the 



intermediate elements (Pd and Ag). Moreover, it is also unable 
to correctly account for the lighter elements (Sr - Zr), where a 
charged particle freeze-out is needed. 

Model II 

The second type of model, which is connected to the low-mass 
SN of a c ollapsing O-Ne-Mg c ore are the two-dimensional (2D) 
models of Wa naio et al.l (1201 ll) . Neutrino interaction is included 
in the model, and the explosion is obtained self-consistently 
without any free parameters. The yields are calculated using 
post-processing networks, in which the output quantities from 
the supernova explosion, such as temperature, density, and Y e 
are applied as inputs without free parameters in contrast to the 
parametrised HEW models. Typical input values are entropy, 
S ~ 10 - 20£fl/baryon (much lower than in the HEW predic- 
tions), and Y e ~ 0.4 to 0.5. The scenario allows neutron-captures 
to take place in the neutron-rich clumps o f matter, which are 
convectively transported to the outer layers. Wan aio et al.l (1201 ll) 
found that there is little production of elements heavier than 
Zr in these conditions but suggest a possible reduction in the 
minimal electron fraction T e ,min below their original value of 
^e.org = 0.4 (because of the limitation in the resolution and the 
two-dimensionality of the model). The imp act that low e r Y e . m ; n 
values have on the yields is explored in Wan aio et al.l ([20 11). 
To test this, as in their work, the additional amount of neutron- 
capture elements created in the neutron-rich clumps with an ar- 
tificially reduced T e .min value is added to their original yields 



Yz(Ye,min) — 



Yz(Y e ,or g )M + ZjYziflAMi 
M + E;AM; 



where AM,- is the relative mass of the elements to be ejected 
at an Y e lower than T ei01g = 0.4, and this mass is set to 2 • 1O~ 5 M . 
The subscript i runs over Y e from T e or g down to F e ,min in steps 
of AY e = 0.005. M is the total mass of the ejecta calculated for 
the higher (original) T e ,org, and this mass is 1.1364O~ 2 M . Yz is 
the predicted yield (mole fraction) of the element with atomic 
number Z. 
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Fig. 26. HD 122 196 (dwarf star) compared to different yields cal- 
culated as a function of T e ,min- The highest electron fractions 
(0.395 and 0.35) are seen to have too few free neutrons to create 
the heavier elements (Ag - Eu). 
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Fig. 27. O-Ne-Mg SN model yields with 0.05 step sized decreasing F e ,min starting from ~ 0.3 to 0.15 compared to observation- 
ally derived abundances of dwarf stars (blue, top) and giants (red, bottom). The abundances were normalised to Sr. These model 
predictions fit the pattern of the r-poor star (HD122563) better than that of the r-rich star (CS 31082-001). 
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Yields calculated with lower i^mm were necessary in order 
to obtain considerable amounts of Pd. 

Wi th electron fractions as small as Y e m ; n = 0.2. IWanaio et al.l 
(1201 lh found that the Ba and Eu abundances in HD 122563 could 
be correctly reproduced. These calculated yields are compared to 
the observationally derived abundances for one metal-poor giant 
in Fig. [26] in which we first constrain the F e ,min values. For sim- 
plicity, calculated yields have been plotted in steps of 0.05. From 
this figure, it is clear that the predictions calculated with an elec- 
tron fraction of 0.395 and 0.35 fail to produce sufficiently large 
amounts of any element heavier than zirconium. Hence, these 
models are not be considered any further. 

The usual seven heavy elements are seen to be created at dif- 
ferent values of T e ,min, where the heavier ones require a lower 
T e ,min because they need more neutrons to be produced. For the 
Sr to Zr abundances, the model predictions calculated with a 
T e ,min in the range from 0.35 to 0.395 show a good agreement 
(i.e. agree within the abundance uncertainty ~ 0.15 - 0.3dex) to 
several of the dwarf stars, Pd and Ag seem to start being pro- 
duced in the proper amounts starting from values of 0.3 down 
to 0.2 (in Te.min)- Similar element-F e m i n relations are found for 
the giants, though the y e m i n values seem to be shifted slightly 
towards higher values (see Fig. |271 and l28li . Bins of the required 
Je.min values needed to create these seven heavy elements are for 
the sake of clarity shown in Fig. [28] 



large amounts of the heavy elements (Ba and Eu) and does not 
in general support the conditions needed for a main r-process. 

7. 1 . Discussion of yield predictions 

For the seven elements here scrutinised, both models provide 
satisfactory agreement with the derived abundances within their 
associated uncertainties in more than 60% of the stars. However, 
neither the HEW predictions nor the O-Ne-Mg SNe model alone 
can provide an agreement with these seven abundances applying 
only one set of input parameters. Since very different entropies 
or electron fractions are needed, two different processes seem 
necessary. The faint O-Ne-Mg SNe could very well be the for- 
mation site of these elements in stars such as HD 122563 and 
HD88609, i.e. stars that are relatively speaking enhanced only 
in the lighter elements or generally speaking r-poor. However, 
these supernovae are not the formation site of abundance pat- 
terns of stars such as CS 31082-001, which is also enhanced 
in main r-process elements such as Eu (r-rich). The entropy or 
neutron-richness in the ejecta from O-Ne-Mg SNe are too low 
to facilitate a main r-process similar to that in the parametrised 
HEW winds explored here. Another possibility for creating the 
r-poor stars would be a HEW with multiple (medium +high) en- 
tropies and low electron fractions (Y e < 0.442) — if possible — 
or a fairly high entropy and an electron fraction of 0.49 cf. Fig. 



Mode II yields 



6.3 
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Giants 
Dwarfs 



Fig. 28. Overview of the yields vs electron fraction. Blue repre- 
sents the dwarfs and red the giants. The atomic number is shown 
on the y-axis. 



This difference in the electron fraction needed to produce Eu 
in dwarf and giant stars, could be due to the different behaviour 
of Eu in dwarfs and giants, or that the giants need large NLTE 
abundance corrections. Since Eu is heavier than Ba, it seems un- 
likely that it would need fewer neutrons (larger T e ,min) to form 
than Ba does. To compare the more extreme cases, we com- 
pare stars with strong r-process enhancements and/or depletions 
(r-rich: CS 31082-001 and HD1 15444, r-poor: HD122563 and 
HD88609) to these model predictions as well as the HEW pre- 
dictions. From Fig. [27] it becomes clear that the O-Ne-Mg core- 
collapse SNe may be the site for stars with weak enhancements 
as seen in both HD122563 and HD88609 (i.e. r-poor), although it 
is clear that the site is insufficiently neutron-rich to produce such 



IWanaio et al ] (120111) showed the effect that higher dimen- 
sional models have on the predicted yields, since ID models 
could not create trans-iron elements below Zr, whereas 2D mod- 
els could (in their Fig. 3). It is therefore important to consider 
3D models before attempting to constrain the r-process site. The 
stellar atmospheres as well as synthesis codes may also need up- 
dates on the input physics, namely improvements in the ID, LTE 
assumptions. With the current model predictions at hand, we 
may need to assume that there are several different sites and/or 
progenitors to explain the diverse abundances patterns we de- 
rive from stellar observations. On the basis of our model com- 
parison, we cannot draw strong conclusions about the weak r- 
process site, but the O-Ne-Mg core-collapse SNe seem promis- 
ing. However, from the abundance pattern of HD 122563 and 
HD88609, the presence of Eu indicates that the ISM must have 
been pre-enriched by main r-process material, or that their gas 
was ejected from an object that allowed both the weak and main 
r-processes to coincide in the same object. For both processes the 
HEW winds are possible sites, since they offer the conditions 
needecPl but it remains unknown whether the high entropies 
are physically feasible. These statements depend very much on 
the efficiency of the mixing in the ISM at metallicities around 
and below about [Fe/H] = -2.5. Therefore, we need to investi- 
gate the abundance patterns of stars in the metallicity range -3.3 
to -5 dex (i.e. below the metallicities of our sample), to deter- 
mine how these patterns that reflect more pristine gases behave 
and compare to the model predictions, before we can address in 
detail the mixing efficiency in the early Galaxy. Unfortunately, 
these abundance patterns cannot be contrasted, since we would 
need very high-quality spectra (S/N > 100 at 3200 A), which 



20 In their current state, the HEW model predictions allowing for a 
large range of parameters, especially the large span in entropies, can 
reproduce several different abundance patterns of all neutron-capture 
elements. Hence, assuming that these physical conditions are viable in 
one single site, t hey explain the patterns with one "continuous broad 
range r-process" (Kratz et al. 2008b; Farouqi et al. 2009; Roederer et al. 
[203). 
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would in turn require observations of very long durations for 
these single faint stars, in order to derive the crucial abundances, 
such as those of Pd and Ag. 

Here as well as in previous studies, it has become evident 
that knowing the precise values of Y e (T e ,min) is essential to ac- 
cu rately predicting the ejected abundance patterns. The work 
of IWinteler et al.l ^20 12) showed that magneto-rotational core- 
collapse SN jets can repro duce the solar abundance pattern for 
120 < A < 210. Moreover. [Arcones & Montesl(l201 ll) argue that 
neutrino-driven winds, either proton- or neutron-rich, stemming 
from core-collapse SNe can create nuclei in the range 65 < A < 
115. Despite the different sites suggested, both studies illustrated 
the importance of knowing Y e . This in turn translates into under- 
standing t he neutrino interactions a nd their effect on the electron 
fraction (Ar cones & Montesll20lTl) . 

From the comparison of our derived abundances to the HEW 
model predictions, we learn that we only weakly detect process 
contributions of the order of 10% or less, but we need to assume 
that there has been more than a 10 - 15% contribution to see 
the effect of the process in our abundance ratios (cf. Table |3)- 
Alternatively, our abundance ratios might actually be sensitive 
down to and below a process contribution of 10%, since some of 
the estimated fractional process contributions might have been 
slightly overestimated. This contribution would change drasti- 
cally depending on the Y e , which is affected by the other free 
input parameters as well as the estimated importance of neutrino 
interactions. 

8. Summary and conclusion 

Summary 

Based on the correlations and anti-correlations of Sr, Y, Zr, Pd, 
Ag, Ba, and Eu, it has become evident with time that the for- 
mation of the heavy elements is not a straightforward process to 
model, and that the previously believed solar-scaled universal r- 
process pattern only provides good estimates for the heaviest el- 
ements, not the elements in the atomic mass range from 38 to 47. 
At low metallicities ([Fe/H] < -2.5), several studies combined 
with this have shown that at least four different neutron-capture 
processes (s- and r-processes with both weak and main com- 
ponents) are needed to explain the observationally determined 
abundances. Starting with the lighter element Sr, which can in 
part be produced via charged particle freeze-out, this element 
might also receive a considerable contribution from a weak s- 
process at higher metallicities. Zirconium has shown similarities 
to both Sr and Ag, which indicates that Zr is created by multiple 
processes, since the formations of Sr and Ag differ. Silver is cre- 
ated by a second/weak r-process. Its formation process clearly 
deviates from the main r-process responsible for forming Eu and 
to some extent Ba at low metallicity. The picture becomes even 
more differentiated when we try to understand the formation of 
Ba. To date, an inexplicably large star-to-star scatter is found for 
the Ba abundances both under LTE and NLTE assumptions. No 
single process can explain all of these results despite the possi- 
bility that some scatter may be caused by the process occurring 
at various sites, in the case of the extremely metal-poor stars. 
Two processes might be needed to fully explain the formation of 
Ba. 

From the comparison to model predictions, we see that de- 
spite the different physics and parameter space investigated, 
the 2-3 dex star-to-star scatter in the stellar abundances cannot 
be explained by e.g. NLTE corrections, stellar parameter influ- 
ence or sample biases, which confirms the need for at least two 



neutron-capture processes yielding heavy elements at very low 
metallicities (below -2.5 dex in [Fe/H]). Some scatter can arise 
from the different amount each site produces, and it seems neces- 
sary to have a combination of various sites to explain the individ- 
ual abundance patterns that the different (r-poor vs r-rich) stars 
show. One possible explanation could be that massive super- 
novae facilitate high entropy winds, which create some amount 
of intermediate elements (in the atomic number interval range 
40 - 50) via a weak r-process, combined with yields of the 
heavi est main r-process elements. According to IWooslev et al.l 
d!994l) . different entropies can be found within one supernova, 
giving rise to various entropies and processes in the exploding 
winds. Another site contributing to the weak r-process elements 
Pd and Ag could be the O-Ne-Mg SNe, which are pred i cted t o 
be very common. However, according to Wanai o et al. 
this kind of SN cannot facilitate a main r-process, hence cannot 
produce Eu, and it seems u nlikely to be the dom inant r-process 
site, which was indicated in llshimaru et al.l(l2004l) . However, 3D 
models might change this picture, though it seems unlikely that 
these supernovae would ever reach the proper conditions to cre- 
ate a full main r-process. Nevertheless, we still need several sites 
and r-processes to explain the abundances of Ag - Eu. 

Conclusion 

We have found that in our observed sample of stars both dwarfs 
and giants show on average the same correlation/anti-correlation 
at all metallicities, thus we feel confident that the correlation 
trends combined with the large star-to-star scatter confirm the 
presence of two different r-processes. A second/weak r-process 
creating Zr - Ag (generally elements in the atomic number range 
40 - 50), and a main r-process producing the very heavy ele- 
ments. The second r-process seems to work in intermediate en- 
tropy, and neutron number density environments, and its path 
possibly lies closer to stability than that of the main r-process. 
However, on the basis of previous studies and the current state 
of the model predictions, we cannot disregard the possibility that 
our suggested second different r-process is in fact a lower end of 
a continuous broad-range main r-process. Many of the physical 
parameters differ between the weak and the main r-process, by 
many orders of magnitude, and we therefore need stronger con- 
straints on what suffices and/or is necessary to define an individ- 
ual process. 

In addition, it seems that several sites are needed to explain 
the diverse stellar abundance pattern coming from r-poor and r- 
rich stars. Possible formation sites are the high-entropy winds 
of SNe and neutron-rich ejecta of electron-capture (O-Ne-Mg) 
SNe. The yields from these objects will be mixed in the ISM, 
which makes tracing the original site a challenging task; further- 
more, these objects are unlikely to be the only sites, and we still 
do not know their frequency or mixing ratios. 
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Appendix A: Line lists 

Here we give the details of our line lists as well as our adopted 
solar abundances and additional useful information for all the 
heavy elements studied. Table I A. il provides the atomic informa- 



Table A.l. Line parameters for the resonance 5s-5p in Agl. 
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-0.936 


107 


Sl/2 


2 Pl/2 


1-1 


3382.885 


-0.635 


109 


Sl/2 


2 Pl/2 


0-1 


3382.894 


-0.936 


109 


2 Sl/2 


2 Pl/2 


1-0 


3382.886 


-0.936 


109 


2 Sl/2 


2 Pl/2 


1-1 


3382.887 


-0.635 


107 


2 Sl/2 


2 P 3 /2 


0-1 


3280.684 


-0.624 


107 


2 Sl/2 


2 ?3/2 


1-1 


3280.678 


-0.925 


107 


2 Sl/2 


2 P3/2 


1-2 


3280.678 


-0.226 


109 


2 Sl/2 


2 P 3 /2 


0-1 


3280.686 


-0.624 


109 


2 Sl/2 


2 ?3/2 


1-1 


3280.679 


-0.925 


109 


"Sl/2 


2 ?3/2 


1-2 


3280.680 


-0.226 



tion for both of the silver lines. These calculations were made 
without any assumption about the natural isotopic ratio. If indi- 
vidual isotopic Ag abundances are needed, these log gf values 
should be applied instead of those listed in Table Q] 

Below we list the solar abund ances that we used. These were 
adopted from lAsplund et all (I2009I) . 

Table A.2. Element and adopted solar abundances. 



Element 


loge 


Sr 


2.87 


Y 


2.21 


Zr 


2.58 


Pd 


1.57 


Ag 


0.94 


Ba 


2.18 


Eu 


0.52 



We provide the atomic data used for the heavy elements 
in ou r line list. The values are taken from VALD dKupka F.I 
2000). The molecular inf ormation in ou r line list was taken from 
Kurucz's home page and Ku ruczl d 1993h . 



Table A.3. Atomic data for the strontium to europium: 
Wavelength, excitation potential, and logg/. The "F indicates 
that th e value is the total l ogg/, which for Ba was s plit accord- 
ing to iMc Willi am (1998) and for Eu according to llvans et al.l 
(200^). 



SrII 


X 






[eV] 


[dex] 


3464.45 


3.04 


0.49 


4077.71 


0.00 


0.17 


4161.79 


2.94 


-0.50 


4215.52 


0.00 


-0.14 


YII 


X 


logg/ 




[eV] 


[dex] 


3549.01 


0.13 


-0.28 


3600.74 


0.18 


0.28 


3628.70 


0.13 


-0.71 


3774.34 


0.13 


0.21 


3788.70 


0.10 


-0.07 


3950.36 


0.10 


-0.49 


4398.01 


0.13 


-1.00 


4854.87 


0.99 


-0.38 


4883.69 


1.08 


0.07 


5087.42 


1.08 


-0.17 


5200.42 


0.99 


-0.57 


Zrll 


X 


logg/ 




[eV] 


[dex] 


3356.09 


0.09 


-0.51 


3499.57 


0.41 


-0.81 


3551.96 


0.09 


-0.31 


3573.06 


0.32 


-1.04 


3607.38 


1.24 


-0.64 


3714.79 


0.53 


-0.93 


4050.33 


0.71 


-1.00 


4161.21 


0.71 


-0.72 


4208.98 


0.71 


-0.46 


4317.32 


0.71 


-1.38 


5112.28 


1.66 


-0.59 


Pdl 


X 


logg/ 




[eV] 


[dex] 


3404.58 


0.814 


0.320 


Ball 


X 


logg/ 




[eV] 


[dex] 


4554.03 


0.00 


0.17' 


4934.08 


0.00 


-0.15 r 


5853.67 


0.60 


-1.01 r 


EuII 


X 


logg/ 




[eV] 


[dex] 


4129.73 


0.00 


0.22' 


4205.04 


0.00 


0.21 7 " 


6645.06 


1.38 


0.12 r 



Appendix B: Stellar parameters 

Below, we list the parameters needed to determine the tempera- 
ture and gravity. Top: Giants, bottom: Dwarfs. The superscripts 
a,b, and c indicate the following: {a) Stars with T e ff and log 
g derived from excitation potential and ionisation balance. 
Stars with a special r-process pattern — either r-poor or r-rich. 
(t) Stellar parameters were altered in accordance with (a) owing 
to uncertainties in colour, dereddening, and parallax. 
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V 


K 


n 
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g 
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O HQ 
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A 
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1 OA 

1.U6 
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U.U1 
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536U 
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2. /6 


O 1 1 

-2.1 1 


l.z 
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O OO 
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l.UU 


n no 

u.yy 
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445U 
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U.5U 
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z.l 


T Tl~"\ 1 1 C A A A b 
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8.96 


6.61 


3.55 


1 11 

1.12 


0.01 




4785 


1.50 c 


-3.00 


2.1 


T Tl~"\ 1 TTCiClfc 

HD1ZZD63 


6.20 


3.73 


3.76 




0.025 




4560" 


0.90" 


-2.81 


1.8 


TT1~"\1 OOriCA 

HD 122*00 


7.25 


5.90 


3.30 




0.083 




4700 


1.51 


-1.45 


1.2 


T Tl~"\ 1 liCIO O 

HD12o23o 


o aa 

7.66 


5.34 


3.81 


0.95 


0.04 




4900 


1.80 


-1.92 


1.5 


T Tl~"\ 1 1 A C O *7 

HDliODo/ 


9.11 


6.668 


1 An 
1.40 


1 /I A 

1.44 


n nn 
0.09 




4700" 


i ncc 

1.05 


O 1 A 

-3.16 


1.7 


T Tl~"\ 1 00^70 

HD12oz/9 


7.97 


7.07 


c nA 

5.96 




n i n 
0.10 




conna 

5200 


o ona 
2.20 


O O /I 

-2.34 


1.3 


HD165195 


O 1 A 

7.34 


A 1 A 

4.14 


2.20 




n 1 nc 

0.195 




a onnc 
4200 


n nnc 
0.90 


o i n 
-2.10 


z.l 


T Tl~"\ 1 S~ S 1 S" 1 h 

HD16616F 


8.12 


5.34 


3.25 


1.19 


0.13 




5250" 


2. 15 c 


-1.25 


1.9 


1 1 U 1 / jjUj 


/.18 


J.UO 


A 1 S 
O.lo 


n ca 

U.5D 


u.uo 




c. 1 nn 
5 1UU 
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z. /u 


1 OS 
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1 

l.Z 


1 1 l_y 1 0U4 / 


Q 14 
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1 TO 
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U.U4 
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4/UU 
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1 0/1 
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Z.U 


JririU j lj+ UUUU 


1 C CO 


1 j.ZU 






n nc i 




conn 
52UU 


o Arv* 
2.4U 


CO 

-z.5y 


1.0 


rlnU44Z-lzj4 


1 o o 1 
12. yi 


O OA 






n 1 11 
U. 1 jj 




/icon 
45oU 
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U.o5 


o oo 
-2.32 


1 
l.O 


tie i o i o 1 o 


1 C Q/1 
15.y4 


1 ^ CO 

1 j.oy 






n nn 




a Ann 
40UU 
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l .U5 
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O.Zl 


1 A 


Star 


V 


K 


7T 


o-(/r) 


E(B-V) 


Mass 


T 


g 


[Fe/H] 




BD+092190 


11.15 


9.91 


1.04 


2.79 


0.0281 


0.8 


6450" 


4.00 


-2.60 


1.5 


BD- 133442 


10.29 


9.02 






0.044 


0.8 


6450 


4.20" 


-2.56 


1.5 


CD-3018140 


9.95 


8.66 
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0.75 
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Notes. w For the giants we assume M = 1M Q 
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Appendix C: Abundances 

Table IC.ll provides an overview of all the abundances and as- 
sociated uncertainties determined for our sample's dwarf stars. 
Table IC.2I gives the abundances and uncertainties for the giant 
stars in our sample. 
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Table C.l. Stellar abundances of Fe, Sr, Y, Zr, Pd, Ag, Ba, and Eu for dwarfs. The ' <' indicates that the abundance is an upper 
limit. 



Star 


[Fe/H] 


[Sr/Fe] 


[Y/Fe] 


[Zr/Fe] 


[Pd/Fe] 


[Ag/Fe] 


[Ba/Fe] 


[Eu/Fe] 


BD+09 2190 


-2.60 


- 


-0.28 ±0.17 


-0.02 ±0.14 


0.72 


±0.27 


- 


- 


- 


BD-13 3442 


-2.56 


0.21 ±0.14 


-0.02 ±0.13 


0.44 ±0.18 




_ 


- 


- 


- 


CD-30 18140 


-1.92 


0.15 ±0.11 


0.1 ±0.12 


0.47 ±0.14 




_ 


- 


-0.10 ±0.25 


- 


CD-33 3337 


-1.55 


0.22 ±0.11 


0.01 ±0.13 


0.27 ±0.17 


0.19 


±0.27 


0.27 ± 0.25 


0.18 ±0.14 


- 


CD-45 3283 


-0.99 


-0.15 ±0.30 


0.03 ±0.15 


0.14 ±0.17 


0.52 


±0.19 


0.38 ± 0.25 


0.32 ±0.15 


0.78 ±0.17 


CD-57 1633 


-1.01 


0.00 ±0.10 


-0.23 ±0.12 


- 


0.17 


±0.23 


0.20 ± 0.26 


0.13 ±0.14 


0.55 ± 0.22 


HD3567 


-1.33 


-0.1 ±0.20 


-0.18 ±0.19 


0.27 ±0.17 


0.29 


±0.20 


0.53 ± 0.25 


0.26 ±0.17 


0.70 ±0.18 


HD 19445 


-2.13 


0.13 ±0.11 


-0.1 ±0.12 


0.29 ±0.15 


0.00 


±0.27 


- 


-0.02 ±0.14 


0.37 ± 0.23 


HD22879 


-0.95 


0.33 ±0.14 


-0.06 ±0.12 


0.19 ±0.17 


0.17 


±0.21 


0.00 ± 0.26 


-0.18 ±0.18 


0.45 ± 0.2 


HD25704 


-1.12 


0.30 ±0.11 


-0.05 ±0.12 


-0.02 ±0.17 


0.07 


±0.21 


0.09 ± 0.26 


0.21 ±0.15 


0.48 ±0.17 


HD63077 


-1.05 


0.36 ± 0.25 


0.06 ±0.13 


0.06 ±0.17 


0.00 


±0.22 


0.07 ± 0.26 


0.22 ±0.15 


0.43 ±0.17 


HD63598 


-0.990 


0.41 ±0.14 


0.09 ±0.12 


0.20 ± 0.22 


0.09 


±0.19 


-0.07 ± 0.26 


0.19 ±0.15 


0.65 ± 0.27 


HD76932 


-0.97 


<0.27 


-0.07 ±0.13 


0.19 ±0.17 


0.20 


±0.20 


0.25 ± 0.25 


0.30 ±0.15 


0.41 ±0.16 


HD103723 


-0.85 


0.04 ±0.13 


-0.27 ±0.15 


0.05 ±0.17 


0.29 


±0.19 


0.21 ±0.25 


0.17 ±0.15 


0.43 ±0.16 


HD 105004 


-0.84 


0.10 ±0.10 


-0.19 ±0.12 


0.02 ±0.17 


0.12 


±0.19 


0.01 ± 0.27 


0.17 ±0.18 


0.26 ± 0.20 


HD106038 


-1.48 


0.56 ±0.17 


0.54 ±0.14 


0.68 ±0.17 


0.32 


±0.21 


0.14 ±0.25 


0.76 ±0.17 


0.45 ±0.16 


HD1 13679 


-0.63 


<0.44 


0.08 ±0.12 


0.22 ±0.17 


0.14 


±0.19 


0.08 ± 0.25 


0.25 ±0.15 


0.08 ± 0.21 


HD 11 1980 


-1.32 


0.45 ± 0.22 


0.23 ±0.13 


0.44 ±0.17 


0.17 


±0.21 


0.10 ±0.25 


0.48 ± 0.15 


<0.5 


HD 11 6064 


-2.17 


- 


0.00 ±0.13 


0.33 ±0.14 


0.37 


±0.21 


- 


-0.36 ±0.17 


- 


HD 120559 


-1.31 


0.22 ±0.12 


0.25 ±0.15 


0.34 ±0.17 


0.44 


±0.21 


0.48 ± 0.25 


0.22 ±0.17 


0.71 ±0.17 


HD121004 


-0.73 


0.40 ±0.12 


0.16 ±0.12 


0.51 ±0.17 


0.32 


±0.19 


0.16 ±0.25 


0.36 ±0.18 


0.50 ±0.19 


HD122196 


-1.81 


0.03 ±0.14 


-0.15 ±0.12 


0.20 ±0.16 


0.02 


±0.27 


<0.22 


0.02 ± 0.16 


0.22 ±0.18 


HD 126681 


-1.28 


0.20 ±0.19 


0.34 ±0.12 


0.56 ±0.17 


0.51 


±0.21 


0.28 ± 0.25 


0.55 ±0.15 


0.47 ±0.17 


HD 132475 


-1.52 


0.34 ±0.11 


0.17 ±0.15 


0.46 ±0.17 


0.34 


±0.19 


0.20 ± 0.26 


0.40 ±0.15 


0.43 ±0.16 


HD140283 


-2.58 


-0.27 ±0.10 


-0.48 ±0.12 


-0.20 ±0.14 






- 


< -0.62 


- 


HD160617 


-1.83 


0.04 ±0.14 


-0.03 ±0.12 


0.19 ±0.21 


0.42 


±0.21 


<0.35 


0.41 ±0.16 


- 


HD166913 


-1.93 


0.47 ±0.15 


0.39 ±0.14 


0.65 ±0.18 


0.42 


±0.23 


<0.63 


0.62 ±0.14 


0.61 ±0.16 


HD 175 179 


-0.72 


<1.28 


<0.95 


0.17 ±0.17 


0.12 


±0.20 


0.09 ± 0.26 


0.45 ±0.14 


0.20 ±0.16 


HD188510 


-1.58 


-0.04 ±0.14 


-0.16 ±0.13 


0.18 ±0.17 


0.32 


±0.19 


0.19 ±0.26 


0.13 ±0.14 


0.40 ±0.18 


HD189558 


-1.18 


-0.70 ±0.14 


0.15 ±0.13 


0.44 ±0.17 


0.30 


±0.21 


0.16 ±0.25 


0.32 ±0.17 


0.39 ±0.16 


HD195633 


-0.71 


<1.11 


<0.66 


-0.14 ±0.17 


-0.13 ±0.20 


-0.03 ± 0.28 


0.10 ±0.16 


<0.1 


HD205650 


-1.19 


-0.02 ± 0.25 


0.05 ±0.13 


0.19 ±0.17 


0.22 


±0.21 


0.14 ±0.25 


0.20 ±0.14 


0.45 ±0.16 


HD213657 


-2.01 


0.04 ±0.11 


-0.05 ±0.12 


0.37 ±0.15 


0.22 


±0.27 


<0.53 


-0.03 ±0.15 




nlJZ70700 


1 zLS 
- 1 .H-o 


U.Uj ± U.l / 


n no + n i a 


v.zj ± u. i / 


0.32 


±0.19 






n ^a -a- n 1 fx 


G 01-039 


-2.27 


0.16 ±0.10 


-0.07 ±0.12 


0.29 ±0.14 












G 05-040 


-0.93 


<1.29 


<1.09 


0.37 ±0.17 


0.16 


±0.21 


0.17 ±0.26 


0.35 ±0.14 


<0.5 


G 20-024 


-1.90 


0.22 ±0.11 


0.17 ±0.12 


0.56 ±0.15 


0.77 


±0.27 




0.32 ±0.14 




G 64-012 


-3.10 


-0.05 ±0.10 


0.02 ±0.13 










-0.35 ±0.16 




G 64-037 


-3.16 


-0.06 ±0.10 


0.03 ±0.13 


0.52 ±0.14 












G 88-032 


-2.53 




-0.14 ±0.13 


0.25 ±0.14 












G 88-040 


-0.89 


0.04 ±0.10 


-0.27 ±0.15 


0.02 ±0.17 


0.12 


±0.19 


0.03 ± 0.26 


0.05 ±0.14 


0.32 ±0.16 


G183-011 


-2.12 




-0.24 ±0.12 


0.18 ±0.17 
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Table C.2. Stellar abundances of Fe, Sr, Y, Zr, Pd, Ag, Ba, and Eu for giants. The ' <' indicates that the abundance is an upper limit. 



Star 


[Fe/H] 


[Sr/Fe] 


[Y/Fe] 


[Zr/Fe] 


[Pd/Fe] 


[Ag/Fe] 


[Ba/Fe] 


[Eu/Fe] 


BD-01 2916 


-1.99 


0.11 ±0.18 


0.03 ± 0.20 


0.26 ± 


0.16 


0.39 ± 


0.20 


0.16 ±0.22 


0.36 ±0.15 


0.60 ± 


0.18 


BD+42 621 


-2.48 


-0.18 ±0.15 


-0.32 ±0.19 


0.12 ± 


0.17 






-0.40 ± 0.32 


-0.55 ±0.14 


0.54 ± 


0.16 


BD+8 2856 


-2.09 


-0.01 ±0.12 


0.04 ±0.15 


0.05 ± 


0.14 


0.12 ± 


0.21 


<0.73 


0.26 ±0.14 


0.41 ± 


0.16 


BD+30 2611 


-1.20 


-0.09 ± 0.26 


-0.35 ±0.18 


0.28 ± 


0.22 


-0.01 ±0.21 


-0.50 ± 0.28 


0.28 ±0.14 


0.52 ± 


0.19 


BD+54 1323 


-1.64 


0.05 ±0.12 


-0.08 ±0.18 


0.01 ± 


0.16 


-0.33 ± 0.20 


<-0.21 


0.35 ±0.14 


0.18 ± 


0.16 


CS22890-024 


-2.77 


-0.06 ±0.18 


-0.25 ±0.16 


0.64±0.16 






- 


-0.30 ±0.16 


<0.9 


CS295 12-073 


-2.67 


0.28 ±0.14 


0.09 ±0.15 


-0.07 ±0.16 






- 


- 


0.08 ± 


0.22 


CS303 12-059 


-3.06 


0.09 ±0.15 


-0.32 ±0.13 


0.43 ± 


0.37 






- 


0.07 ±0.16 






CS30312-100 


-2.62 


-0.34 ± 0.20 


-0.79 ±0.12 


-0.17 ±0.17 






- 


- 


0.12 ± 


0.21 


CS3 1082-001 


-2.81 


0.66 ±0.11 


0.82±0.16 


0.77 ± 


0.26 


1.29 ± 


0.19 


1.19 ±0.27 


1.43 ±0.14 


1.53 ± 


0.31 


HD74462 


-1.48 


0.06 ±0.10 


0.38 ±0.16 


0.54 ± 


0.26 


0.13 ± 


0.18 


-0.33 ± 0.25 


0.37 ±0.14 


0.52 ± 


0.17 


HD83212 


-1.25 


-0.04 ±0.17 


0.21 ± 0.22 


0.27 ± 


0.22 


-0.20 i 


= 0.19 


-0.53 ± 0.28 


0.32 ±0.14 


0.27 ± 


0.17 


HD88609 


-2.87 


0.04 ±0.11 


-0.09 ± 0.23 


0.19 ± 


0.14 


<0.07 


<0.25 


-0.93 ±0.14 


-0.50 = 


t0.16 


HD108317 


-2.11 


-0.05 ±0.13 


-0.22 ±0.12 


0.07 ± 


0.14 


-0.08 i 


= 0.20 


0.15 ±0.32 


0.34 ±0.17 


0.34 ± 


0.18 


HD1 10184 


-2.40 


-0.05 ± 0.22 


0.15 ±0.14 


0.47 ± 


0.18 


0.22 ± 


0.23 


0.15 ±0.22 


-0.06 ±0.14 


0.22 ± 


0.16 


HD 11 5444 


-3.00 


-0.08 ±0.16 


-0.12 ±0.14 


0.20 ± 


0.19 


0.47 ± 


0.20 


<0.35 


0.35 ±0.14 


0.79 ± 


0.16 


HD122563 


-2.81 


0.04 ±0.13 


-0.10 ±0.16 


0.07 ± 


0.18 


0.12 ± 


0.20 


0.10 ±0.25 


-0.85 ±0.14 


-0.51 ±0.16 


HD122956 


-1.45 


-0.04 ±0.14 


0.00 ± 0.28 


-0.01 ±0.25 


0.32 ± 


0.20 


-0.28 ± 0.22 


0.34 ±0.14 


0.24 ± 


0.16 


HD 12623 8 


-1.92 


-0.09 ±0.14 


-0.27 ±0.12 


0.06 ± 


0.14 


0.24 ± 


0.19 


-0.01 ± 0.28 


0.16 ±0.14 


0.19 ± 


0.16 


HD126587 


-3.16 


-0.01 ±0.10 


-0.20 ±0.13 


0.19 ± 


0.19 


0.42 ± 


0.26 


0.40 ± 0.28 


0.08 ±0.14 


0.31 ± 


0.17 


HD128279 


-2.34 


-0.36 ±0.11 


-0.78 ±0.13 


-0.35 = 


t0.14 


-0.16 = 


= 0.20 


-0.24 ± 0.22 


-0.43 ±0.15 


-0.28 = 


t0.16 


un i sir i (iC 


-2.10 


—0.19 ± U.lO 


—0.01 ± U.2U 


0.09 ± 


0.20 


-0.32 i 


= 0.20 


-0.6U ± U.23 


U.5S ± U.14 


0.89 ± 


0.17 


HD166161 


-1.25 


0.29 ±0.14 


0.29 ±0.18 


0.26 ± 


0.17 


0.00 ± 


0.19 


0.05 ± 0.25 


0.55 ±0.15 


0.07 ± 


0.16 


HD175305 


-1.38 


0.01 ±0.19 


0.12 ±0.16 


0.29 ± 


0.18 


0.04 ± 


0.18 


0.04 ± 0.28 


0.35 ±0.15 


0.56 ± 


0.16 


HD 186478 


-2.42 


0.08 ±0.11 


0.01 ±0.12 


0.31 ± 


0.14 


0.26 ± 


0.18 


0.23 ± 0.25 


0.25 ±0.14 


0.50 ± 


0.16 


HD204543 


-1.84 


-0.07 ±0.12 


-0.14 ±0.15 


0.07 ± 


0.18 


0.04 ± 


0.18 


0.04 ± 0.25 


0.23 ±0.15 


0.10 ± 


0.17 


HE 0315+0000 


-2.59 


0.13 ±0.12 


0.08 ±0.15 


0.30 ± 


0.16 










0.70 ± 


0.21 


HE 0442-1234 


-2.32 


-0.18 ±0.17 


-0.24 ±0.16 


-0.06 ±0.19 


0.12 ± 


0.30 






0.28 ± 


0.14 


HE 1219-0312 


-3.21 


0.01 ±0.13 


-0.19 ±0.21 


0.19 ± 


0.14 








0.70 ± 0.20 







